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Abstract 

We review the hadro-production data presently available on open charm and 
beauty absolute production cross-sections, collected by experiments at CERN, 
DESY and Fermilab. The published charm production cross-section values are 
updated, in particular for the "time evolution" of the branching ratios. These 
measurements are compared to LO pQCD calculations, as a function of the colli- 
sion energy, using recent parameterisations of the parton distribution functions. 
We then estimate, including nuclear effects on the parton densities, the charm and 
beauty production cross-sections relevant for measurements at SPS and RHIC 
energies, in proton-proton, proton-nucleus and nucleus-nucleus collisions. The 
calculations are also compared with measurements of single D and B kinemati- 
cal distributions, and DD pair correlations. We finish with two brief comments, 
concerning the importance of beauty production as a feed-down source of J/V' 
production, and open charm measurements performed using leptonic decays. 
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1 Introduction 



In the context of the study of heavy ion colhsions, heavy flavour production is be- 
coming increasingly more interesting, as the energies available for particle production 
increase from fixed target (SPS) to collider (RHIC, LHC) experiments. Charm is the 
heaviest flavour which can be produced in nucleus-nucleus collisions at SPS energies, 
where experiments with high intensity beams and a dimuon trigger made detailed 
studies of the production and suppression of charmonium states (J/'?/', 'ip'). Besides 
being the natural reference for charmonia studies (same initial state), open charm 
significantly contributes to the yield of dimuons measured in the mass range be- 
tween the and the ^/ip resonances, through the simultaneous semi-muonic decays 
of a pair of D mesons. Enhanced production of continuum dileptons in this mass 
range, as seen by several experiments at the CERN SPS, could be a signal of thermal 
dimuon production from a Quark-Gluon Plasma phase, emphasising the importance 
of understanding the charm "background" . 

At the higher energies of RHIC and LHC, also beauty production comes into play, 
not only as a direct probe of the properties of the very early phase of the collision 
system but also as a source oi ^/ip mesons. Indeed, the very important study of J/?/' 
production (and suppression or enhancement) at collider energies requires a good 
understanding of the fraction of ^/ip^s coming from the feed-down of B decays. If 
this feed-down source is not well understood, no proper physics interpretation of the 
data will be possible. 

Therefore, it is very important to have a good evaluation of the yield of charm 
and beauty production as a function of the collision energy. This is necessary for the 
understanding of the SPS data and to make reliable estimates of the yields expected 
at the collider experiments. To ensure their usefulness, these evaluations should be 
made using a well known and easily available calculation procedure, such as the one 
provided by the event generator Pythia. However, it is essential to ensure that the 
generator properly describes existing data, collected over the last years, essentially in 
fixed target experiments with proton and pion beams. 

Since the calculations critically depend on the distribution functions of the quarks 
and gluons present in the interacting hadrons, we must see how the results vary if 
we use different sets of parton distribution functions (PDFs). Finally, since we want 
to use these evaluations in the context of nuclear collisions, it is also important 
to correctly estimate the effects of the nuclear modifications of the PDFs on the 
production of heavy flavours. However, the degree of nuclear shadowing at low values 
of X is not well known for the valence and sea quark distributions, and there is no 
direct information at all on the nuclear effects on the gluon distribution function, at 
any value of x. This makes the predictions of heavy flavour production particularly 
uncertain for nuclear collisions, especially at the higher energies (lower x values) 
available at the colliders, where gluon fusion dominates the production cross-sections. 

In the following section we will discuss some basic concepts relevant to the topic 
of heavy flavour hadro-production, with emphasis on the importance of the parton 
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distribution functions, including nuclear effects. In Section 01 we introduce the ex- 
periments which provided the data included in this report. The data are critically 
reviewed and compiled in Section |^ for open charm production and in Section El for 
the beauty case. These data are then compared, in Section IHl to the LO pQCD 
calculations provided by the event generator Pythia. In this section we also estimate 
the cross-sections for experiments made at SPS and RHIC energies, and we compare 
the calculations with existing data on single meson kinematical distributions and on 
pair correlations. We finish with brief comments on the relevance of the beauty feed- 
down to the studies oi ^/ip suppression at RHIC and higher energies, and on some 
(indirect) measurements of charm production yields using leptonic decays. 

We should clarify that several issues related to heavy-flavour production are not 
addressed in this review. In particular, we only addressed hadro-production data; 
see, for instance, Ref. P (and references therein) for information related to photo- 
production or e^e~ experiments. We have also not covered in here any effects related 
to "heavy-ion physics", such as charm flow, gluon saturation, recombination, dead- 
cone effect, etc; see Refs. 0121, for instance, for details on such issues. And it is clear 
that we only addressed open heavy flavour production; see, for instance, Refs. [HI El 
for details on quarkonia production, in elementary and nuclear collisions. 

2 Heavy flavour production in pQCD 

In this section we will briefly review some basic issues related to the physics of heavy- 
flavour production in hadronic collisions. For a more detailed introduction, see e.g. 
Ref. p. 

At leading order (LO), the only processes which can lead to heavy flavour pro- 
duction are quark- ant iquark annihilation and gluon fusion, as illustrated in Fig. ^ 
In general, pQCD processes can be factorised into three different parts: the non- 
perturbative initial conditions, describing the state before the collision takes place, 
the hard process itself, perturbatively calculable, and the subsequent step of hadro- 
nisation (also non-perturbative) . 




Figure 1: Heavy flavour production mechanisms at leading order. 

The first part is mainly determined by the interacting quarks' and gluons' frac- 
tional momenta, a;i 2 = Pparton/Pi,25 where 1 and 2 stand for the projectile and target 
nucleons. It also depends on the squared energy- momentum transfer between the two 
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partons, Q^, and on the intrinsic transverse momenta, k^, which the partons carry 
inside the projectile (proton, pion or nucleus). The distributions of the fractional 
momenta of the various quarks inside protons and pions were studied in deep inelas- 
tic scattering (DIS) experiments, among others, and were parameterised by several 
groups. They will be discussed in the next section. The "strength" of the pro- 
cess is given by the partonic cross-section, aij, which depends, in particular, on the 
available energy. At a given energy, in LO calculations, the partonic cross-section is 
determined by the mass of the heavy quark, mq, and by the strong coupling constant, 
as, evaluated at the scale Q^, 

ffj^ is a dimensionless scaling function which determines the energy dependence of 
the charm or beauty production cross- sect ion, and which depends on the ratio mq/s, 
where s is the squared partonic centre of mass energy, s = X1X2S [S]. The indices 
represent the interacting partons {ij = qq or gg). 

The cross-section to produce a heavy quark pair in a proton-proton collision, 
0"^?^, is then obtained by convoluting the perturbatively calculated partonic cross- 
section with the (non-perturbative) parton distribution functions, /p, of the interact- 
ing hadrons, 

^QQ = E / ■ ■ /f (xi, Q') ■ ff{x2, Q') ■ a.,{s) . (2) 

If the protons are inside nuclei, their partons have modified distribution functions. 
In Eq. (0) the parton distribution functions should then be represented by f^{x, Q^) 
instead of f^{x,Q^), where A represents the mass number of the nucleus. The im- 
plications of this effect on the total cross-sections will be studied in Sections 12.21 
and 

In the last step, the hadronisation, the heavy quark pair fragments into hadrons, 
including the neutral and charged D or B mesons. We will discuss this step in 
Section 1231 

2.1 Parton distribution functions 

Figure |21 shows the fractional momentum distributions, in protons, of the valence 
quarks, sea quarks and gluons. The valence quarks usually have relatively high mo- 
mentum fractions while the sea quarks and the gluons are mostly found at low x 
values. As first observed by NA51 [HI and then studied in more detail by E866 [7j, 
the u and d distributions are not identical; at x ~ 0.2 there are almost twice as many 
d than u in the proton sea. The right panel of Fig. |21 shows how the probabilities 
of finding certain partons at given momentum fractions change when we probe them 
in reactions of different values. The higher the energy- momentum transfer, Q^, 
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Figure 3: Various sets of proton LO (left) and NLO (right) PDFs. 



the higher the probabihty that low-x partons participate in the particle production 
process. 

Later in this report we will calculate the dependence of cc and bb production, 
using different PDF sets. In view of this work, we have upgraded the PDFLIB 
package jH] to include the recent CTEQ6 [0] and MRST 2001 (TU] parameterisations. 

Figure El shows different sets of LO and NLO proton PDFs. While the four LO sets 
show a similar behaviour, the NLO parameterisations show a more significant spread, 
in particular for the gluons. The CTEQ6M parameterisation differs from other sets 
mainly in the following two aspects: an extended function is used to fit the data 
points, including correlated systematic errors, and new measurements are included. 
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with improved precision and expanded {x, Q"^) ranges. Among the new data sets, the 
inclusive jet cross-section measurements of the D0 experiment at Fermilab, giving 
access to the x range 0.01 < x < 0.5, are particularly important, since they have 
a big impact on the CTEQ6M gluon distribution functions. In Ref. jH] the CTEQ 
Collaboration gives a detailed description of their new method, which mainly uses 
the MS scheme, and outlines the differences with respect to the MRST 2001 sets. 

At the bottom of Fig. IHl we roughly indicate the ranges probed by fixed target 
(SPS, FNAL, DESY) and collider (RHIC, LHC) experiments, for charm and beauty 
production. They were evaluated using the expression x = M/^/s ■ exp(y*), with 
masses of 1.5 and 5 GeV/c^ for charm and beauty, respectively. We set y* = 0, where 
y* is the rapidity in the centre of mass frame, and varied the energy within the ranges 
of the available experimental measurements. 




Figure 4: Three different sets of NLO pion PDFS (left) and comparison of pion and 
proton LO PDFs (right). 

In this report we will consider measurements of open charm and beauty production 
cross-sections from data obtained with proton and pion beams. In the left panel of 
Fig. m we show three pion PDFs, calculated at NLO. It should be noted that all 
available pion PDF sets are more than 10 years old. These three pion PDFs are 
significantly different from each other. On the right panel of this figure we compare 
LO parton distribution functions in a pion and in a proton. We can see that the 
valence quark distributions are peaked at x ~ 0.2 in protons and x ~ 0.45 in pions, 
where the gluons are much harder. 

It is clear that the use of different PDF sets will change the calculated total 
production cross-section of heavy flavour production (see Eq. HJ. It also influences 
the relative importance of gluon fusion and qq annihilation. Figure El shows the 
relative contribution of gluon fusion to the total cc production cross-section, as a 
function of ^/s, as calculated by the Monte Carlo event generator Pythia , for pp 
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Figure 5: Relative contribution of gluon fusion to the total cc production cross- 
section, as a function of -^/i, in pp (left and middle) and 7r~p (right) collisions. 
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Figure 6: Same as previous figure but for bb production, and up to higher energies. 



and 7r~p collisions. The remainder of the total cross-section is due to qq annihilation, 
since at LO there are only these two processes. More details on the calculations will 
be given later. Note that for the 7r~p collisions we used the CTEQ6L (2002) and 
MRST c-g (2001) sets to describe the parton distributions inside the target proton, 
but other sets of proton PDFs give similar results. Most measurements of the cc 
cross-section were made in the range 200 < i?iab < 920 GeV, or ^/s = 20-40 GeV, 
where the contribution from gluon fusion is around 80 % in pp and around 70 % in 
7r~p collisions. 

In Fig. iniwe show the y/s dependence of the relative contribution of gluon fusion 
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to the total bb cross-section. The higher values obtained with the CTEQ6M and 
GRV98 HO PDFs, at the lowest energies, result from the harder gluon distributions 
of these sets. Comparing Figs. El and IHl we see that in the energy range of the 
fixed target experiments the production of the two heavy flavours is dominated by 
different mechanisms: while qq annihilation is responsible for only ~ 20 % of the total 
cc production cross-section, it is the dominant process in bb production. At higher 
energies, both charm and beauty production are dominated by gluon fusion. 



2.2 Nuclear effects in p-A and A- A collisions 



The parton distribution functions in the proton, /P(a;,(5^), are essentially extracted 
from the structure functions (Fi, F2 and F3) measured in deep inelastic scattering 
experiments. These experiments are performed with various nuclear targets, and 
indicate that the distributions of partons inside bound protons are different from 
those in hydrogen. These nuclear effects are expressed as the ratio of the PDFs 
observed in a nucleus with respect to those of a "free" proton. 



with i representing the valence quarks, the sea quarks, or the gluons. 
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Figure 7: Nuclear modifications of the PDFs, for the Pb nucleus, according to the 
EKS 98 [12] weight functions. 

Figure [7| shows this "nuclear weight function", for the Pb nucleus, as a function of 
X, according to the EKS 98 ^21 parameterisation. The curves are shown for two values 
oi . A detailed discussion of the dependence of the nuclear effects on Q^, among 
many other related issues, can be found in Refs. ^21 El- The interested reader will 
find in Ref. ^31 a recent and detailed review of nuclear parton distribution functions. 
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In different regions of x, the nuclear effects are traditionally referred to by the 
following expressions: 

• shadowing: at low x, where R^{x,Q^) < 1 

• anti- shadowing: at medium x, where R'^{x,Q^) > 1 

• EMC effect: at relatively high x, where R^^x^Q"^) < 1 

• Fermi motion: at the highest x, where R'^i^x^Q'^) > 1 




Figure 8: Changes induced on the cc and bb cross-sections by the nuclear modifica- 
tions of the PDFs, at mid-rapidity. 

The impact of these nuclear effects on heavy flavour production can be seen in 
Fig.lHl which shows the ratio between the heavy flavour (cc and bb) production cross- 
sections calculated for pp collisions taking into account that those protons are inside 
Pb nuclei (p-Pb, Pb-Pb) and the same cross-sections calculated for pp collisions in 
the vacuum. According to the EKS 98 parameterisation, and in what concerns charm 
production, the experiments carried out at the SPS and FNAL energies are in the 
anti- shadowing regime. Therefore, inside a heavy nucleus the parton distributions 
are expected to be harder than in a "free" proton, leading to higher production 
cross-sections in p-A and A-A collisions with respect to the linear scaling from pp 
collisions. For instance, at ^/s ~ 30-40 GeV, the cc production cross-sections in 
p-Pb collisions should be 10% higher than in the absence of gluon anti-shadowing. 
The RHIC experiments, at mid-rapidity, are just in the x range where the nuclear 
effects on charm production change from the anti-shadowing to the shadowing region. 
Therefore, measurements in the central detectors should not be very sensitive to 
nuclear effects on the parton distribution functions. However, this is no longer the 
case for the detectors placed away from mid-rapidity. At ^/s = 200 GeV, and for 
charm production, we have the same gluon anti-shadowing (~ 15%) at y = —2.0 as 
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we have at the SPS {y/s = 27.4 GeV) at mid-rapidity. In the case of d-Au collisions, 
there is a significant difference between the expected nuclear effects on the PDFs 
in the "North" and "South" muon arms of Phenix, for instance. In the "forward 
hemisphere" (with respect to the d beam), at y = +2.0, there is ~20% shadowing 
effects on charm production, instead of the ~ 15 % anti- shadowing expected on the 
"backward" side. It is crucial, hence, to keep the PHENIX charm and charmonia d-Au 
analyses independent for each of the three covered rapidity ranges. In what concerns 
beauty, the nuclear effects are expected to influence the production cross-sections 
in opposite ways when going from measurements done in the energy range of fixed- 
target experiments (EMC region) to those done at RHIC (anti-shadowing region), 
with the nuclear cross-sections changing from suppressed to enhanced, with respect 
to the linear extrapolation from pp collisions. At the LHC energies, we will certainly 
be in the shadowing region, both for charm and for beauty. Charm production in the 
Pb-Pb collision system, for instance, is expected to be suppressed by around 40 % 
with respect to a linear extrapolation of nucleon-nucleon collisions. 

It should be noted, however, that there are no measurements today which can con- 
strain the nuclear gluon distribution function. The nuclear gluon densities provided 
by the EKS98 parameterisation are only indirectly constrained, through the scale evo- 
lution of F^{x, Q"^) and through momentum conservation. Accurate measurements of 
open charm production in proton-nucleus collisions, using several nuclear targets and 
over a broad range of energies, would be crucial to significantly reduce the present 
uncertainties on the nuclear gluon densities [THj . 



2.3 Fragmentation 

In the hadronisation step, the outgoing heavy quarks fragment into hadrons. The 
energy carried by the formed hadron with respect to the quark's energy, z = Eh/Eq, 
is distributed according to the fragmentation function, Dq{z), measured in e~^e~ 
reactions and assumed to be the same in hadronic collisions. We should note that 
this definition of z is not unique; in theoretical studies it is more common to use the 
lightcone fraction, where E is replaced by E + ph, with the longitudinal momentum 
defined by the quark's direction. Light quark (u,d,s) fragmentation is usually param- 
eterised as D^{z) oc z~^{l — z)" while heavy quarks only experience a relatively small 
deceleration when combining with a (slow) light quark. In 1983, Peterson et al. p3] 
proposed the following heavy quark fragmentation function 

which peaks at 2; ~ 1 — ^/eg. In principle, eg is fixed by the light and heavy quark 
masses, eg = (mg/mg)^, but in practice it is a free parameter, usually taken to be 
0.06 for charm and 0.006 for beauty. 

Alternatively, the Lund string fragmentation scheme ^H] leads to the expression 

Dq{z)<^ ^\'^^^J -expi-bml/z) , (5) 

z Q 
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which is sensitive to the of the produced hadron, through the transverse mass. It 
imphes a harder fragmentation function for heavier quarks, through the exphcit mass 
dependence. The term z-''"'Q was introduced by Bowler ^Tj in 1981, to improve the 
agreement with the e'^e~ B meson data available at that time. 




Figure 9: Heavy flavour fragmentation functions according to the Peterson ^3] and 
Lund string ^H] expressions. All curves are normalised to unity. 

Figure ini shows both fragmentation functions, for charm and beauty hadrons. The 
Lund model curve was drawn with a = 0.3, b = 0.58 GeV~^, rric = 1.5 GeV/c^ and 
TTT-b = 4.8 GeV/c^ (values used by the Pythia event generator [H]), and with an 
exponential spectrum, dN/dpT oc exp(— 5 ■ p^). 



3 Heavy flavour experiments 

Most of the experiments considered in this report were designed to study the proper- 
ties of charm and/or beauty hadrons. To select events with charm or beauty particles, 
the experiments used high-resolution detectors in the target region, to observe pri- 
mary and secondary vertices, signaling the decay of the heavy flavoured hadrons. 
Three types of vertex detectors have been commonly used: bubble chambers, emul- 
sions and silicon tracking telescopes. The spacial resolution of the silicon detectors 
is not as good as that of the other systems, but they can be operated at much 
higher interaction rates, a crucial feature when looking for rare processes. Most of 
the experiments had particle tracking devices and a muon spectrometer. In addition, 
the experiments which measured the charm (or beauty) hadrons in hadronic decay 
channels also had electromagnetic and/or hadronic calorimeters, and particle identifi- 
cation detectors, such as Cerenkov counters, to distinguish pions, kaons and protons. 
Furthermore, many of them, in particular those which studied beauty production, im- 
plemented triggers to enrich their collected event sample with charm and/or beauty 
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events. With only one exception (E789), the fixed target experiments could detect 
the charm or beauty hadrons in the full forward hemisphere. Some experiments had 
active targets, where the vertex detector itself was used as target. 

The heavy-fiavour hadro-production experiments we consider in this report are 
summarised in Table Q and described in some detail in the next pages, roughly in 
chronological order. 



Flavour 


Experiment 


charm 
beauty 
both 


NAU, NA16, NA27, E743, NA32, WA75, WA82, E769, E791, E706 
NAIO, WA78, UAl, E706/E672, E771 
E653, E789, WA92, CDF, HERA-B 



Table 1: Heavy-fiavour hadro-production experiments considered in this report. 



3.1 NA16, NA27 and E743: the LEBC Experiments 

The purpose of NA16 was the "study of hadronic production and properties of new 
particles with a lifetime 10"^^ < r < 10"^° s, using LEBC-EHS". It used the high- 
resolution hydrogen bubble chamber "LEBC" (Lexan bubble chamber) and a proto- 
type version of the European Hybrid Spectrometer, "EHS" fH]- LEBC was a rapid 
cycling liquid hydrogen bubble chamber, with a fiducial volume of 12 x 5 x 2.5 cm'^ 
photographed by two cameras. It served both as a liquid-hydrogen target and as a 
high-resolution vertex detector. The direct observation of the production and de- 
cay vertices is one of the key features of this apparatus. The decay products were 
analysed downstream in the EHS spectrometer, which could detect photons but had 
very limited particle identification capabilities. The acceptance for D mesons covered 
the positive xp range and was independent of the observed decay mode. The data 
samples were collected in the late seventies, at the CERN SPS, with 360 GeV proton 
and 7i~ beams |T^ . 

NA27 was built "to measure accurately the lifetime of the D*^, D^, F''^, A^ charm 
particles and to study their hadronic production and decay properties" (note that 
the F''^ is now named D^). It used the final version of the EHS spectrometer |2Uj . 
composed of three parts and extending over more than 40 metres. Immediately 
downstream of LEBC were placed two wire chambers, complemented by two small 
drift chambers, for track reconstruction and triggering purposes. The trigger simply 
required more than two hits in each of the wire chambers. Each of the two other parts 
had a magnet and three large drift chambers, leading to a momentum measurement 
with a relative resolution better than 1 % up to 250 GeV/c. Electron and photon 
detection were provided by two lead-glass electromagnetic calorimeters. A hadron 
calorimeter was also available. Charged particle tracking and identification were 
essentially performed by a 40 m^ drift chamber, ISIS, through up to 320 dE/dx 
measurements, complemented by two Cerenkov detector systems and a transition 
radiation detector. In the 400 GeV proton run, a total of 98 neutral and 119 charged 
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D mesons were found [21]. Previously, NA27 had a vr" run at 360 GeV, with lower 
statistics [22] • The data samples were taken in the early eighties, at the CERN SPS. 

Some of the CERN ISR experiments studied charm production indirectly, by trig- 
gering on single electrons. The deduced cross-sections were between ten and hundred 
times higher than those observed in the LEBC experiments. Since the 800 GeV pro- 
ton beam of Fermilab provided collisions with an energy half way between the SPS 
and the ISR data, an experiment was proposed, E743, to investigate the discrepancy 
between the previous measurements. The bubble chamber LEBC was transported to 
Fermilab and complemented with a multi-particle magnetic spectrometer, MPS, that 
had Cerenkov and transition radiation detectors for particle identification, besides 
proportional wire chambers for tracking. The interaction trigger was provided by 
two proportional wire chambers placed just downstream of LEBC. The experiment 
collected data in 1985. Like NA16 and NA27, also E743 identified the charm mesons 
in topological decays, by observing the charge of the decaying particle and a given 
number of charged final state particles |2S]. The new measurement (10 neutral and 46 
charged D mesons) agreed with the results of the previous SPS LEBC experiments. 

3.2 NAIO 

NAlO was designed to perform a "high resolution study of the inclusive production of 
massive muon pairs by intense pion beams" and took data in the early eighties, at the 
SPS [21]. The muon spectrometer, separated from the target region by a 5 m long 
carbon muon filter, was composed of eight multi-wire proportional chambers with 
three tracking planes each, and four trigger hodoscopes made of plastic scintillator 
slabs, separated in two telescopes by an air core toroidal magnet. The last trigger 
hodoscope was protected by a 1.2 m iron wall, placed after the tracking chambers, 
to ensure a clean dimuon trigger without deteriorating the reconstruction of the 
muon trajectories. The highly selective dimuon trigger, optimised for masses above 
3 GeV/c^, allowed to run at a beam intensity of ~ 1.5 x 10^ pions/burst. The study of 
BB production in 7r~-W interactions was based on the selection of events with three 
high pt muons in the final state, coming from the semi-muonic decays of both B 
mesons and from the semi-muonic decay of one of the D mesons. Beauty production 
cross-sections were given [22] for incident beam energies of 140, 194 and 286 GeV, 
the largest statistics (~ 14 signal events) corresponding to the highest energy. 

3.3 WA78, WA75 

WA78 was proposed to "search for the hadropro duct ion of BB pairs" and took place 
in the early eighties, at the SPS [22]. It followed a similar strategy as NAIO, looking 
at three muons in the final state or at like-sign muon pairs. In addition to the muon 
spectrometer, consisting of drift and multi-wire proportional chambers surrounding 
a 1.5 T superconducting dipole magnet, WA78 had extensive calorimetry. Because 
of the large mass difference between beauty and charm mesons, muons produced in 
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the B decay have larger and are accompanied by more energetic neutrinos than 
those produced by charm decays. The trigger and event selection procedures were, 
therefore, designed to select events with at least two high-px muons and large missing 
energy. The final event samples, collected with a 320 GeV vr" beam incident on an 
U target, contained both tri-muon events and like-sign muon pairs fI7\ . 

The same muon spectrometer had previously been used by WA75 j2H]; which 
collected a few hundred charm events (and one BB event). The online event se- 
lection required at least one high-px muon, to enhance the fraction of events with 
semi-leptonic decays of heavy fiavour particles. The primary and secondary vertices 
were located in the emulsion target, within the volume indicated by the tracks re- 
constructed with the silicon microstrip planes. A total of 339 events were observed 
with the identified muon among the tracks of the decay vertices, presumably due to 
charmed particle semi-leptonic decays. 

3.4 NAll and NA32: the ACCMOR Experiments 

The ACCMOR Collaboration started their charm physics program in 1980, as NAll, 
taking 7r~-Be data at 120, 175 and 200 GeV j^H], at the SPS. Their large acceptance 
forward magnetic spectrometer included two magnets, four sets of drift chambers 
(a total of 48 planes), a complex system of multi-wire proportional chambers, five 
Cerenkov counters and Pb-scintillator electromagnetic calorimeters. The data sam- 
ples were collected with a single electron trigger. In a second stage, a vertex tele- 
scope of high resolution silicon microstrip detectors was added to the setup, helping 
to reduce the background levels. The statistics of the charm event sample, however, 
remained small [301 1^ • 

In 1984, the NA32 experiment took over, to "investigate charm production in 
hadron interactions using high-resolution silicon detectors". The ACCMOR spec- 
trometer [221 complemented with a finely segmented active silicon target, made 
of 14 planes of 20 /im pitch silicon microstrip detectors, 280 fim thick, preceeded by 
a silicon beam telescope and followed by two silicon multiplicity counters. Forward 
going particles were tracked in 7 planes of 20 yum pitch silicon microstrip detectors, 
before entering the ACCMOR spectrometer. The data samples were taken with an 
interaction trigger, defined using signals from the active target. Most of the statistics 
came from runs with a 200 GeV hadron beam [23| consisting of vr" and K~, separately 
identified by means of two threshold Cerenkov beam counters. Much fewer events 
were also collected with a 200 GeV proton beam, allowing the experiment to publish 
cross-sections for three different beam particles. 

The second stage of the NA32 experiment was performed in 1985/86 with a 2.5 mm 
copper target, placed in vacuum, an improved vertex detector, including CCDs of 
22 X 22 yum^ pixels, and a two-level trigger to select Ac and Ds decays [21]. In these 
runs, the 230 GeV hadron beam contained 96 % vr" and 4 % K~ mesons. 
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3.5 E653 



The E653 experiment was designed to "study charm and beauty using hadronic 
production in a hybrid emulsion spectrometer" . It was the first experiment designed 
to measure both charm and beauty hadrons, and took data in 1985 (with protons) 
and 1987 (with pions), at Fermilab. It used a 1.47 cm long emulsion target with 
A = 26.6 as the average nuclear mass number. This allowed to measure the primary 
vertex and at least one decay vertex still within the emulsion volume. The trigger, 
which was optimised to select semi-muonic decays of charm particles, required an 
interaction in the target and a high muon candidate. Tracking started with a 18- 
plane silicon microstrip vertex detector and continued with a magnetic spectrometer 
composed of a dipole magnet and 55 drift chamber planes. After a 5 m long steel 
absorber was the muon spectrometer, made of 12 drift chamber planes on each side 
of an iron toroidal magnet. The primary vertex was located visually, and the muon 
trajectory was compared to the tracks in the emulsion in order to find a good match. 

Cross-sections for neutral and charged D meson production were obtained with 
800 GeV protons jSHl and 600 GeV pions jHZj, essentially on the basis of semi-muonic 
decays, with a small contribution from purely hadronic decay channels. For the 
beauty study [2H1; the 600 GeV pion data sample was carefully re-analysed to look for 
additional vertices. As before, the emulsion analysis procedure selected events with 
muonic secondary vertices of large muon transverse momentum, before proceeding 
with the search for other decay vertices. 9 events were found where both the B and 
the B decay vertices could be identified and the whole decay chain reconstructed. 

3.6 E769, E791 

The open geometry TPL ("tagged photon laboratory") spectrometer, at Fermilab, 
after being used by E516 and E691 to measure charm photo-production, was "in- 
herited" by the E769 experiment to study "pion and kaon production of charm and 
charm-strange states", in 1987/88. The E769 detector included an 11-plane silicon 
microstrip vertex detector, 2 analysing magnets, 35 drift chambers, 2 multi-wire pro- 
portional chambers, and 2 segmented threshold Cerenkov counters to identify kaons, 
pions and protons. It also included electromagnetic and hadronic calorimeters, a wall 
of scintillation counters for muon identification, and a high-rate data acquisition sys- 
tem. E769 used a mixed 250 GeV secondary beam of both charges. A much smaller 
data sample was also collected with a 210 GeV negative beam. The composition of 
the negative beam was 93% vr", 5% K~ and 1.5% p, while the positive beam was 
61% 7r~^, 4.4% K"*" and 34% p. The beam-particle identification was provided by 
a differential Cerenkov counter complemented, in the case of the positive beam, by 
a transition radiation detector. Eight proportional wire chambers and two silicon 
microstrip planes were used to track the beam. This allowed the E769 Collabora- 
tion to study D meson production in pion, kaon and proton induced collisions |39j . 
The mixed hadron beam collided on a multifoil target, consisting of 250 /^m thick 
Be, Al and Cu foils, as well as 100 fim thick W foils, interspaced by 1.6 mm. In 
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total, 26 target foils were used, giving a total target thickness of 2 % of an interac- 
tion length. Having different target materials simultaneously in the beam allowed 
the experiment to measure the nuclear dependence of open charm production in pion 
induced collisions. 

E791 was approved to study "hadropro duct ion of heavy flavours at the tagged 
photon laboratory". With respect to the E769 detector, the E791 experiment in- 
creased the number of silicon microstrip planes in the vertex telescope, added a sec- 
ond scintillator wall in the muon identifier, and implemented a faster read-out system. 
The 2 % interaction length target was made of one 0.52 mm thick platinum and four 
1.56 mm thick diamond disks, interspaced by 1.53 cm. This rather large spacing 
ensured that the decay of a charm hadron would occur between the target foils. It 
collected 88 990 neutral D mesons in 1991/92, with a pure 7i~ beam of 500 GeV [10] . 

3.7 UAl 

The UAl detector was built to find the intermediate vector bosons in pp collisions, 
using the SPS in a collider mode. It was basically composed of a cylindrical drift 
chamber and an electromagnetic calorimeter immersed in a dipole magnetic field, 
surrounded by a hadron calorimeter and a 8-layer muon detector. "End-cap" electro- 
magnetic and hadronic calorimeters were installed in the forward directions, giving 
the detector an excellent hermetic coverage. The beauty production cross-section 
measurement |1J was performed in the 1988 and 1989 runs, at ^/s = 630 GeV, when 
the central electromagnetic and forward calorimeters were removed in preparation 
for the installation of new detectors. The muon detection system was improved by 
the addition of iron shielding in the forward region. Muon trigger processors selected 
tracks in the muon chambers pointing back to the interaction region. At high lumi- 
nosity, the muon trigger rate in the forward region was further reduced by requiring 
a jet of transverse energy greater than 10 GeV in coincidence with the muon trigger. 
The search for beauty hadrons was performed in four independent decay channels: 
b — > fiX, b — i> J/ip{^ /i+/i~)Ar, b c(— > /iX) fiX and b — > fiX;h fiX. Cross- 
sections were measured for each of these processes, each channel covering different 
ranges in from which B hadron and b quark cross-sections were inferred. The 
combined cross-section was then extrapolated to full phase space. 

3.8 E672 and E706 

E706 was designed to perform "a comprehensive study of direct photon production in 
hadron induced collisions", at Fermilab. The detector complemented a large accep- 
tance liquid argon calorimeter, containing a finely segmented electromagnetic section 
and a hadronic section, with a charged particle tracking system composed of sili- 
con microstrip detectors, a large aperture dipole magnet, proportional wire chambers 
and straw tube drift chambers. The experiment collected events triggered by high 
transverse momentum showers detected in the electromagnetic calorimeter. This re- 
quirement enhanced the fraction of selected events containing charm by nearly an 
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order of magnitude, compared to a minimum bias trigger. The measurements |42j . 
restricted to charged D mesons, were performed in 1990, using a negative 515 GeV 
beam, primarily composed of pions with a small admixture of kaons (< 5 %), not sep- 
arated. Two 780 fim thick copper targets were followed by two beryllium cylinders, 
3.71 and 1.12 cm long. 

Downstream of the E706 apparatus, 20 m away from the target and protected by 
a steel wall to absorb most of the hadrons, was placed the E672 muon spectrome- 
ter, aimed at "studying hadronic final states produced in association with high-mass 
dimuons" . The E672 muon spectrometer was composed of six proportional wire cham- 
bers, a toroidal magnet and two scintillator hodoscopes, besides iron and concrete 
shielding, to provide a clean dimuon trigger. In 1990, the E706 and E672 Collabora- 
tions joined efforts to study beauty production in 7r~-Be collisions at 515 GeV 
using J/ip mesons coming from secondary vertices to tag the beauty candidates. The 
trigger selected dimuons in the proximity of the J/ip mass; they had a mass resolution 
of 68 MeV and an average vertex resolution of 14 and 350 fim, in the transverse and 
longitudinal coordinates, respectively. 

3.9 E789 

E789 was proposed to "measure the production and decay into two-body modes 
of b-quark mesons and baryons", and took data in 1990/91. E789 upgraded the 
spectrometer previously used by the E605 and E772 experiments by adding a vertex 
telescope, made of 16 planes of 50 fim pitch microstrip silicon detectors and placed 
between 37 and 94 cm downstream of the target, to identify the decays of neutral 
D mesons. Data samples were taken with either beryllium or gold targets. A vertex 
processor selected (on-line) track pairs consistent with decay vertices at least 1.02 mm 
downstream of the target and impact parameters of at least 51 /xm relative to the 
target centre. The spectrometer featured two large magnets. Particles were identified 
by electromagnetic and hadronic calorimeters, scintillation hodoscopes, proportional- 
tube muon detectors and a ring-imaging Cerenkov counter. E789 collected a large 
statistics data sample |13], much bigger than all other charm measurements made 
with proton beams, but only measured neutral D mesons and their acceptance was 
limited to the < xp < 0.08 window. It also measured the nuclear dependence of 
neutral D meson production, comparing data taken on beryllium and gold targets. 

E789 was the first experiment to measure the beauty production cross-section in 
proton-nucleus interactions using a gold target of 50 x 0.2 mm^ area and 3 mm 
thick. The target was placed in vacuum, to ensure that interactions in air would 
not be confused with 6-hadron decays. Within its acceptance, the highly energetic 
(piab ~ 150 GeV) B hadrons had their production and decay vertices separated by 
an average distance of 1.3 cm. E789 triggered on events with a dimuon coming from 
the target region, to look for beauty hadrons through their decay into a. J/ip meson. 
The spectrometer had excellent dimuon mass resolution: 16 MeV at the J/ip mass, 
dominated by multiple scattering in the target. The longitudinal vertex resolution 
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was 700 /iin. J/?/^'s from beauty decays were requested to have their origin more than 
7 mm downstream of the target centre. The impact parameter of each muon, defined 
as the vertical distance between the muon track and the target centre, had to be 
larger than 150 fim. 19 ± 5 events survived these rather strict selection cuts. 

3.10 E771 

The E771 Collaboration upgraded the Fermilab High Intensity Lab spectrometer, 
previously used by E537 and E705, with silicon microstrip detectors, pad chambers 
and resistive plate counters, to "study charm and beauty states as detected by decays 
into muons" jlHl- The beauty production cross-sections were measured from data 
collected during one month in 1991, using the 800 GeV proton beam. The target 
consisted of twelve 2 mm thick Si foils interspaced by 4 mm, giving a total effective 
length of 5.2 % Aint- A silicon microvertex detector was positioned downstream of the 
target for the measurement of primary and secondary vertices. Multiwire proportional 
chambers and drift chambers were used, together with a dipole analysis magnet, to 
determine charged particle trajectories and momenta. The spectrometer finished 
with a muon detector made of three planes of resistive plate counters, embedded 
in steel and concrete shielding. E771 used a dimuon trigger to select two possible 
decay modes: a. J/ip coming from the decay of a B meson, or a muon pair from 
the simultaneous semi-muonic decays of two beauty hadrons. Muons from the semi- 
muonic decay of B mesons could be accepted if within the —0.25 < xp < 0.50 window. 

3.11 WA92, WA82 

The Beatrice Collaboration, WA92, used the Omega Spectrometer, at the CERN 
SPS, to "measure beauty particle lifetimes and hadroproduction cross-sections", but 
also published results on charm production ^[Sj. It took data in 1992 and 1993 with 
a 350 GeV vr" beam (including a 1.2% K~ contamination) incident on Cu and W 
targets. Charged particle tracking started with a series of high granularity (10 to 
50 /im pitch) silicon-microstrip detector planes, organised in a Beam Hodoscope, a 
Decay Detector and a Vertex Detector. Charged particles were then tracked in multi- 
wire proportional chambers placed inside the superconducting Omega dipole magnet, 
with a bending power of 7.2 Tm, and in drift chambers placed downstream of the 
magnet. The setup was completed with an electromagnetic calorimeter, followed by 
a muon identifier made of resistive plate counters protected by hadron absorbers. A 
multi-component trigger was developed to identify events with beauty decays, which 
also slightly enriched the fraction of collected charm events. Due to the high statistics 
collected with the Cu and W targets, WA92 also measured the nuclear dependence of 
neutral and charged D meson production The beauty production cross-section 
measurement was derived from the 7r~-Cu data samples jSO]- 

Part of these detectors had already been used by the WA82 experiment, where the 
Omega Spectrometer was complemented by 23 silicon-microstrip detector planes, with 
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a pitch ranging from 10 to 50 fim, to precisely reconstruct tracks and (secondary) 
vertices jSIj. The fast onhne treatment of the sihcon-microstrip data was used to 
select events with at least one track missing the primary vertex. To study the nuclear 
dependence of charm production [52] , the experiment took data (with a 340 GeV vr" 
beam) using a 2 mm thick target made of two materials (either Si/W or Cu/W), 
placed side by side, transversely with respect to the beam axis. The beam illuminated 
simultaneously both target materials, reducing the systematic uncertainties. 

3.12 CDF 

CDF is a "general purpose" experiment, at Fermilab, which studied pp collisions 
at y/s = 1.8 TeV between 1992 and 1995 ("Run I"), and at 1.96 TeV from year 
2001 onwards ("Run II"). Charged track trajectories are reconstructed in a drift 
chamber and matched to strip clusters in the silicon vertex detectors. These devices 
are immersed in a magnetic field of 1.4 T, generated by a superconducting solenoid. 
The central muon system (outside a hadron calorimeter) consists of eight layers of 
drift chambers, four before and four after a 60 cm thick steel absorber, and detects 
muons with > 1.4 GeV/c in the range \y\ < 1. The measurements of D and 
B meson production are based on rather large data samples, but the results from 
Run I are restricted to a relatively high window, > 5.5-6.0 GeV/c, at mid- 
rapidity, \y\ < 1. In Run I, CDF studied charged B meson production using 
the B+ ^ decay channel. Charged B candidates were selected by combining 

the J/ip mesons with each charged particle track of pt > 1-25 GeV/c (kaons from B 
meson decays have a harder spectrum than most other particles). The dimuon 
and kaon tracks were then constrained to come from a common vertex. This study 
was done using events triggered on two opposite-sign muons in the mass range of the 
J/ip, corresponding to an integrated luminosity of Ci^it = 98 ± 4 pb~^. 

Run II provided much better data for heavy flavour studies, largely because of the 
added ability to trigger on secondary vertices. Results on charm production [SI] were 
obtained from data collected in early 2002. The D mesons were selected by requiring 
two oppositely charged tracks with p^ > 2 GeV/c and px.i + Pt,2 > 5.5 GeV/c, 
with a distance of closest approach to the beam axis between 120 /xm and 1 mm. 
Results on beauty production have also been published [SH], based on data collected 
in 2002, with C-mt = 39.7 ± 2.3 pb~^. Beauty hadrons were measured in the inclusive 
B— >■ J/ijj X channel, using the transverse distance between the J/ip origin and the pp 
collision vertex, within \y\ < 0.6 and down topx = 0, a rather impressive experimental 
achievement: close to = the signal drops fast but the background does not, the 
acceptance changes rapidly, etc. 

3.13 HERA-B 

The HERA-B fixed-target experiment, at DESY, was designed to identify, within a 
large geometrical coverage, the decays of B and J/ip mesons produced on target wires 
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by the halo of the 920 GeV HERA proton beam. Most of the events were collected 
with C and W targets, with only a small fraction (less than 10%) taken with Ti. 
Primary and secondary vertices were reconstructed by the vertex detector system, 
made of double-sided silicon microstrip detectors integrated in the HERA proton 
ring. The main tracking system was placed upstream of a 2.3 T dipole magnet. The 
inner region, near the beam pipe, used microstrip gas chambers, while the outer 
tracker was made of honeycomb drift cells. Particle identification was performed by 
a ring imaging Cerenkov detector, a muon spectrometer composed of four tracking 
stations, and an electromagnetic calorimeter. In the data taking period of 2002/2003, 
164 • 10^ events were collected with a dilepton ^/ip trigger. Decays of J/^/', V^', Xc and 
T mesons were reconstructed, in the and in the e^e~ decay channels. 

A bb production cross-section measurement was published jSH] from data collected 
in year 2000, with dimuon and dielectron triggers, through the study oih ^/ip 
I'^l" decays in the window —0.25 < x^^ < 0.15, on the basis of 1.9^^5 dimuon 
events and 8.6132 dielectron events. Improved results from the 2002/2003 run were 
recently published jSZl; following the same analysis, yielding 46.2l7g and 36.9I78 
events in the muonic and electronic channels, respectively. Results on open charm 
production are expected to be published soon. Preliminary values have already been 
made pubhc [SH]- 



4 Data on open charm production 

In this section we collect and discuss the experimentally measured production cross- 
sections for the charged and neutral D mesons, (t(D"'") -|- cr(D~) and o"(D°) -|- (j(D*^). 
Within the last 30 years various experiments, using different kinds of detectors, have 
collected data on open charm production. In the late seventies, four experiments 
at the ISR pp collider, at CERN, reported results on charm production, mostly 
triggering on single electrons, assumed to come from the semi-electronic decay of a 
D~ or a D*^. We have not included these measurements in our study. In some cases, 
only upper limits or ranges were given for the cross-sections. Moreover, due to lack 
of statistics, sometimes data collected at y/s = 52 and 62 GeV were merged to obtain 
a common result on the production cross-section. Besides, the published values differ 
significantly between the different experiments, as discussed in Ref. |59j . 

There are other early experiments jHSj which studied open charm production but 
which we will not consider here, such as NA18, NA25, E515 and E595. They could 
not separate the different charm hadrons, only giving "associated charm production 
cross-sections". Within their rather large uncertainties, their values are consistent 
with the measurements we have considered. We have also ignored the result obtained 
at ^/s = 630 GeV with a modified UA2 detector [HI] , given its huge uncertainty. 

In Table 121 we summarise the data used in the present study, obtained with proton 
and pion beams, at energies ranging from i^iab = 200 to 920 GeV. Very significant sta- 
tistical samples have been collected by WA92 and E791 (~ 7000 and ~ 90 000 events. 
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Table 2: Experiments measuring the production cross-sections of neutral and charged 
D mesons. The average mass number of the emulsion target of E653 is A = 26.6. 



respectively) with pion beams, while the fewer proton beam experiments collected 
much less data. The ~300 events of E769 (adding neutral and charged D mesons 
and the statistics of four different nuclear targets) constitute the highest statistics 
proton event sample, among the fixed target experiments covering a reasonably large 
phase space window. 

At the much higher energies of the Tevatron pp collider, ^/s = 1.96 TeV, CDF 
collected (within the \y\ < 1 window) 36 804 D*^ mesons, of p^ > 5.5 GeV/c, and 
28 361 D+ mesons, of pt > 6.0 GeV/c. 

In order to properly compare the different measurements to each other, we applied 
certain corrections to some of the published values. In particular, whenever possible, 
we normalised the published cross-sections to the latest branching ratio values, using 
the Particle Data Group 2004 values jHl] • Table El summarises the decay channels of 
D mesons which were investigated in each experiment. It further lists the branching 
ratio used in the original publication, together with the corresponding latest value. 
Figure ITIIl illustrates the "time evolution" of two branching ratios. 

We have also updated the systematic errors of the published values to reflect 
the smaller uncertainties of the most recent branching ratios. Some publications 
have not included on their systematic errors these uncertainties, something we must 
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Table 3: Comparison between the latest branching ratios [02] for the various D me- 
son decay channels and the ones used in the original publications. The number in 
parentheses indicates the year of publication. Since the branching ratio of the decay 
channel DO — > K (892)~/i"'"i/^ is not known, the electronic value was used by E653, 
assuming lepton universality. 



do when comparing D meson production measured in different decay channels. If 
the D mesons were searched for in more than one decay channel, the performed 
corrections were weighted according to the number of observed events in each of 
the decay channels. This procedure was not applied to the data of experiments 
which searched the D mesons in topological decays, where the search is not done 
in a specific decay channel, but rather by detecting a certain number of charged or 
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Figure 10: Time evolution of two relevant branching ratios. 
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Table 4: Published and updated D meson production cross-sections in proton in- 
duced collisions. The NA27 errors represent the combined statistical and systematic 
uncertainties. The HERA-B values are preliminary (see the "Note added in proof" 
at the end of this paper). 



neutral final state particles (NA16, NA27, WA75, E743). The neutral D mesons 
were searched, e.g., in "V2" or "V4" prong decays, corresponding to neutral decays 
with 2 or 4 prongs, while "C3" is a typical topological decay of a charged D meson 
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Experiment 


-E'lab 


o given 


o [/ib] 






[GeV] 


for 


published 


updated 


vr- beam, a{D^) + (t(DO) 


NA16 (84) 


m 


360 


Xf > 


7 7+7.2 
'•'-3.5 


7 7-+7.2 
'•'-3.5 


NA27 (86) 


122] 


360 


xf > 


10.1 ± 2.2 


10.1 ± 2.2 ± 15% 


NAll ('86) 




200 


a]l tt? 


31 ± 7 ± 16 


45 ± 10 ± 23 


NA32 (88) 


ESI 


200 


xf > 


1.15 (3.3l[]i± 0.3) 


4.41°:^ ± 0.4 


NA32 (91) 


El 


230 


xf > 


6.3 ±0.3 ± 1.2 


6.6 ±0.3 ±1.0 


E653 (92) 


EH 


600 


xf > 


22.05 ± 1.37 ±4.82 


18.86 ± 1.17 ±3.80 


E769 (96) 


ES 


210 


xf > 


6.4 ±0.9 ±0.3 


6.8 ±1.0 ±0.3 


E769 (96) 


m 


250 


xf > 


8.2 ±0.7 ±0.5 


8.7 ±0.7 ±0.6 


WA92 (97) 


m 


350 


Xf > 


7.78 ±0.14 ±0.52 


7.83 ±0.14 ±0.48 


E791 (99) 


EOl 


500 


Xf > 


15.4+1-8 

— z.o 


15.6+^-8 








vr beam 


, o\\J ) + o\\J ) 




NA16 (84) 


CHI 


360 


Xf > 


A c:+2.2 
^•0_1.4 


4 Ci+2-2 
^•0_1.4 


NA27 (86) 


122] 


360 


Xf > 


5.7 ±1.5 ±0.4 


5.7 ± 1.5 ± 1.4 


NAll (86) 


EO] 


200 


all Xf 


20 ± 5 ± 10 


25 ± 6 ± 12 


NAll (87) 


El] 


200 


all Xf 


30.0 ± 3.5 ± 12.6 


21.3 ±2.5 ±7.3 


NA32 (88) 


EH] 


200 


Xf > 


1.15 (1.7+[J:^± 0.1) 


1.9lJ5|±0.2 


NA32 (91) 


El 


230 


Xf > 


3.2 ±0.2 ±0.7 


2.7 ±0.2+°:^ 


E653 (92) 


EI] 


600 


Xf > 


8.66 ± 0.46 ± 1.96 


8.05 ± 0.43 ± 1.69 


E769 (96) 


Ei 


210 


Xf > 


1.7 ±0.3 ±0.1 


1.7 ±0.3 ±0.1 


E769 (96) 


Ei 


250 


Xf > 


3.6 ±0.2 ±0.2 


3.6 ±0.2 ±0.3 


WA92 (97) 


HH] 


350 


Xf > 


3.28 ±0.08 ±0.29 


3.24 ±0.08 ±0.28 


E706 (97) 


112] 


515 


Xf > 


11.4±2.7±3.3 


11.3±2.7±3.3 


7r+ beam, aijf) + (t(DO) 


E769 (96) jHi 


250 


Xf > 


5.7 ±0.8 ±0.4 


6.0 ±0.8 ±0.4 








TT"*" beam 


, a(D+)±a(D-) 




E769 (96) 


250 


Xf > 


2.6 ±0.3 ±0.2 


2.6 ±0.3 ±0.3 



Table 5: Published and updated D meson production cross-sections in 'pion induced 
reactions. The two charged D meson values of NAll correspond to different decay 
channels (see Table EJ). 



with 3 prongs. Since topological decays are composed of different decay channels, 
with unknown contributing fractions, we cannot correct them for the evolution in our 
knowledge of the branching ratios. Tables ID and El collect the published and updated 
values of the production cross-sections for fixed-target experiments using proton and 
pion induced collisions, respectively. 

CDF measured in Run II j^l , the single D meson production cross-sections in pp 
collisions at = 1.96 TeV, in \y\ < 1: a(D^) = 13.3±0.2±1.5 /ib for > 5.5 GeV/c 
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and a(D+) = 4.3 ± 0.1 ± 0.7 for pt > 6 GeV/c. 

With the exception of WA75, the experiments which used nuclear targets assumed 
a linear dependence on the mass number of the target nucleus to derive the cross- 
section in pp (or Trp) collisions, (TpA = A ■ CTpp. 

We will now comment on certain specific measurements. 

• NA16 (84) The publication does not mention if the quoted errors include sys- 
tematic uncertainties, due to the branching ratios or to other factors. In any 
case, the quoted errors are very large (and certainly dominated by statistical 
uncertainties) and increasing them would not change the overall comparison 
with the other, more precise, data points. Therefore, we have not updated the 
NA16 values or error bars. This applies to all NA16 measurements, with proton 
or pion beams, producing charged or neutral D mesons. 

• E789 (94) The acceptance of the E789 detector was limited to the ranges 
< xp < 0.08 and < 1-1 GeV/c. It is clear that the extrapolation of 
the measurement to full phase space is model dependent, and a significant un- 
certainty should be added to the published systematic error, before comparing 
it to other measurements or to calculations. 

• NA27 (86) The neutral D mesons were searched in the 2-prong and 4-prong 
decays. By weighing the contributions from these two modes, we deduced an 
uncertainty of 15 % in the cross-section measurement, due to the branching 
ratio. This value was used to estimate the systematic error of the published 
cross-sections, since no other value was explicitly quoted in the publication. 

• NA27 (86) The published systematic error of the D charged cross-section is 
0.4 /ib, a value quoted as being mostly due to uncertainties in the branching 
ratio. However, the uncertainty due to the branching ratio, on its own, would 
already give an error of 1.3 /ib. We presume that there was a misprint in the 
original publication and updated the systematic error to 1.4 /ib. 

• NA32 (88) In a later publication [21] the NA32 Collaboration explains that a 
normalisation error was found and all total cross-sections of Ref. |2S| should be 
upscaled by 15 %. 

• NA32(88), NA32 (91) The neutral D mesons are identified in two decay 
channels. Their relative contribution was estimated according to 

jydecays _ ^2body _|_ ^4body _ £ . _ ^^2body _ ^2body _|_ ^4body _ jj4body'j 

where 2body and 4body stand for the 2- and 4-body decay channels, listed in Ta- 
ble El N'^'^^^^^ is the total number of reconstructed D mesons, C the integrated 
luminosity, the cross-section for D meson production, A the acceptances, 
and B the branching ratios. The relative populations of 2-body and 4-body 
decays are then determined to be 39%/61 % for NA32 (88) and 23%/77% for 
NA32(91). The corrections on the branching ratios were then performed ac- 
cording to these weights. 
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• E653 (92) Several decay channels contributed to the detection of the neutral 
D mesons, but we only corrected the branching ratios of those contributing the 
most, listed in Table El 

While most of the results obtained with proton beams were published for the 
full x-p range, those obtained with pion beams were published for the positive xp 
range, with the exception of NAll, which assumed a symmetric xp distribution for 
the extrapolation. This measurement should be taken with care. NAll was the 
predecessor of NA32 and used a "prototype version" of the experimental apparatus. 
NA32 used significantly upgraded detectors and performed a much more accurate 
measurement, with the same pion beam and at the same collision energy. 



Experiment 


Elah 

[GeV] 


a 
[/ib] 


p beam, full xp 


NA32 (88) 


200 


2 (1.5 ± 0.7 ± 0.1) 


E769 (96) 


250 


11.2± 1.7±0.8 


NA16 (84) 


360 




NA27 (88) 


400 


18.1 ± 1.7 


E743 (88) 


800 


29^^ ± 5 


E653 (91) 


800 


48±6± 11 


HERA-B (05) 


920 


51.7 ±5.8 ±6.6 


IT 


beam. 


xf > 


NAll (86) 


200 


21 ±3±8 


NA32 (88) 


200 


3.8lJ^J±0.3 


E769 (96) 


210 


5.1 ±0.6 ±0.2 


NA32 (91) 


230 


5.6 ±0.2 ±0.7 


E769 (96) 


250 


7.4 ±0.4 ±0.4 


WA92 (97) 


350 


6.64 ±0.10 ±0.33 


WA75 (92) 


350 


7.5±0.4+J:? 


NA16 (84) 


360 


7.3ty 


NA27 (86) 


360 


9.5 ± 1.6 ± 1.2 


E653 (92) 


600 


16.15 ±0.75 ±2.50 


IT 


beam. 


xf > 


E769 (96) 


250 


5.2 ±0.5 ±0.3 



Table 6: cc production cross-sections in p-A and vr-A collisions. See the text for 
remarks on the NA32 proton value and on the WA75 value. The HERA-B value is 
preliminary (see the "Note added in proof" at the end of this paper). 

In Table El we summarise the derived values for the total cc production cross- 
sections. The p-A values are given for full phase space while the tt-A results are 
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given for the positive xp hemisphere. To obtain the total cc production cross-section, 
besides adding the measured neutral and charged D meson values, we must take into 
account the production of other charmed hadrons (Ds, Ac and other charmed baryons 
and charmonia states). Only NA32 and E769 also measured the Dg and Ac hadro- 
production cross-sections [311123 IM], but with very poor statistics. Assuming that 
the fragmentation fractions are universal, we can use e"'"e~ or 7p data. World averages 
of results obtained by CLEO, ARGUS, HI, ZEUS and the four LEP Collaborations 
show that neutral and charged D meson production represent 78 ± 3 % of the total 
charm production cross-section This value is well reproduced by the Pythia 

Monte Carlo event generator. Note that the D meson cross-sections compiled in 
Tables 0] and El refer to the sum of the particle and anti-particle values. 

The proton measurement of NA32 was not listed in the previous tables because 
only the total D meson cross-section was published not separating the charged 
and neutral states. The quoted number was scaled up by the factor 1.15 mentioned 
above. The factor 2 explicitly shown in the table converts the published positive 
Xp value to full xp, for consistency with the other proton beam measurements. The 
rather big statistical uncertainty is due to the fact that only 9 events were observed. 
Also WA75 published [2H1 a charm production cross-section without separating the 
neutral and charged contributions, from 339 pion-emulsion events observed in the 
window —0.5 < xp < 0.8. The published vrp value, 23.1 ± 1.313;° A^b/nucleon, for full 
Xp, was derived assuming an A" scaling with a = 0.87. The value included in TablelHl 
was recalculated assuming a = 1.0. It was also divided by 2 so that it roughly 
corresponds to the a^p > window, for consistency with the other measurements. 
The published systematic errors of these two measurements, by NA32 and WA75, 
are underestimated, since they ignore the uncertainty on the branching ratios. Not 
enough information exists in the original publications to calculate a proper error or 
to update the values using the most recent branching ratios. 

5 Data on open beauty production 

The available measurements on beauty production were collected over the last 15 
years. Besides the fixed target experiments, with energies within 200 < i?iab < 
920 GeV, UAl and CDF measured the beauty cross-section at the much higher ener- 
gies of pp colliders, -y/i = 630 GeV and 1.8-1.96 TeV, respectively. Since fixed-target 
experiments barely have the minimum energy required for beauty production, the 
cross-sections are very small, only a few nb, and very selective triggers are needed 
in order to observe even a few beauty events. In most cases this was achieved by 
triggering on high pt single muons or dimuons. Looking at inclusive muonic decays, 
as was done by the older experiments, has the advantage of exploring rather large 
branching ratios, ~10%. The new experiments explored the beauty decay into J/i(j 
by measuring the fraction of J/tp^s with a minimum offset with respect to the inter- 
action point, which has the advantage of having a well understood reference process 
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for normalisation purposes and is much more robust with respect to backgrounds: B 
decays are the only source of displaced i /ip mesons, while there are several sources of 
single muons. These experiments, hence, complemented a muon spectrometer with 
a vertex detector to reconstruct the secondary vertices, profiting from the B mesons' 
long lifetimes. All experiments, except one, reported cross-sections for a global mix- 
ture of all different beauty hadrons, B"*", B~, B°, B°, B^, B°, b-baryons, etc, with a 
priori unknown relative fractions. The exception is the single B"*" meson cross-section 
obtained by CDF in Run I, from a measurement of charged B mesons. 



Experiment 


Beam E^y, 
[GeV] 


iarget 


Phase space window 
(Pt in GeV/c) 


BB events 


NAIO 




7r~ 


286 


W 


> 0; selection cuts 


14 


WA78 


EH 




320 


U 


xp > 0; selection cuts 


73 


E653 






600 


emulsion 


xf > -0.3 




9±3 


E672/E706 


I43j 


TT^ 


515 


Be 


xf > 




8 ±3.3 


WA92 


[■50i 


7r~ 


350 


Cu 


xf > -0.2 




26 


E789 


m 


P 


800 


Au 


0<xJ/'^<0.1, 


4/^<2 


19 ±5 


E771 


m 


P 


800 


Si 


-0.25<x^<0 


.50, p^>6 


15 


HERA-B 


[57J 


P 


920 


C,Ti,W 


-0.35 < Xp/"^ 


< 0.15 


83 ± 12 


UAl 




PP 


v^ = 


= 630 GeV 


\y\ < 1-5, pt 


> 6 


2859 


CDF Run I 


m 


PP 


v^ = 


= 1.8 TeV 


\y\ < 1-0, p^ 


> 6 


387 ±32 


CDF Run II 




PP 


v^ = 


1.96 TeV 


\y\ < 0.6 




38078 



Table 7: Beauty hadro-production measurements. E653 used an emulsion target (of 
average mass number 26.6). In Run I, CDF only measured B''^ mesons. 



Table [7| shows the available measurements of beauty production, the phase space 
covered and the number of identified BB events. Rather than updating the older 
measurements with the most recent knowledge on branching ratios, an issue of minor 
importance when compared to the other uncertainties involved and the large statisti- 
cal errors, we prefer to emphasise that many experiments relied on theoretical models 
to determine the published cross-sections. This is necessary because of the limited 
phase space coverage of the detectors. Often very different models were used, giving 
different extrapolation factors to obtain the total (full phase space) cross-sections. 
Only the E653, WA92, E771 and HERA-B experiments cover around 90 % of full 
phase space and, thus, are less sensitive to theoretical assumptions. 

The first beauty hadro-production cross-section measurement, by NAlO, was pub- 
lished in 1988 j2ni, before QCD calculations were developed in detail, so that they 
used a rather simple production model. NAIO also provided another cross-section 
value, lower by 30 %, obtained by using different xf and pt distributions for the 47r 
extrapolation. The second beauty hadro-production measurement, made by WA78, 
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is in good agreement with the NAIO result, if the same production model is used |27j . 
However, when WA78 uses a LO QCD calculation the two results differ significantly, 
even though the measurements were performed at almost the same energy. 

During the past 20 years, the theoretical understanding of beauty production has 
significantly evolved and the experiments profited from this progress when analysing 
their data. Since they were performed over a long time span, each experiment used 
a different calculation (Pythia Monte Carlo generator, NLO calculations with several 
settings, etc), and it is highly non-trivial to "normalise" all the available measure- 
ments to one common production model. Such an update would require knowing in 
detail all the kinematical cuts used and having the measured values before extrapo- 
lating to full phase space, something not published by all experiments (or we would 
need to start from the published full phase space cross-sections and undo the extrapo- 
lations made in the original analyses, if we could have access to the codes used at that 
moment). It is always better when the experiments publish their results in a form 
which closely reflects the measurements made, in terms of phase space window and 
particles measured. Extrapolations to full phase space, as well as "deconvolutions" 
of the experimental data to the bare heavy quark level, will be biased by theoretical 
prejudice 



Decay channel 


Experiment 


BR used [%] 


BR (PDG04) [%] 


b^fiX 


UAl 


10.2 ± 10% 


inn 


10.95+°:! 


b J/ipX, 


UAl 
E672/706 


1.12 ±0.18 
1.3 ±0.2 ±0.2 


EH 

m 


1.16 ±0.10 




E789 


1.30 ±0.17 








UAl 


6.9 ±0.9 




5.88 ±0.10 




E672/E706 
E789 


5.90 ±0.15 ±0.19 
5.97 ±0.25 


[H] 




B+ ^ J/ipK 


CDF I 


(10.9 ± 1.0) X 10-2 


m 


(10.0 ±0.4) X 10-2 


{b DfiXf 


NAIO 
WA78 


11.32 
11.62 


[IS] 


(11.95 ±0.56)2 


D ^ /iX 


NAIO 


11 




9.2 ±0.8 




WA78 


10.4 


m 





Table 8: Decay channels of the identified BB events and corresponding branching 
ratios, as used in the original publications and in the PDG 2004 tables. In the 
E672/E706 case, although the combined uncertainty of the branching ratios is 22%, 
only 13 % was taken into account in the calculation of the systematic error bar. 



Table IHl collects the various decay modes in which the B mesons were looked for, 
with the corresponding branching ratios, as used by the different experiments and 
as given in the PDG 2004 tables. HERA-B and CDF II used the PDG 2004 values, 
while E771 used the default values of Pythia. E653 is not listed since it reconstructed 
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± lidov; oUdL-c wiiiLivjvv 










(px ill GeV/c) 




[nbl 

I J 


NAIO (88) 










WA78 (89) 


m 





— 


3.6 ± 0.4 ± 1.1 


E653 (93) 


m 




— 


33 ±11 ±6 


E672/706 (95) 


m 


xp > 


47 ± 19 ± 14 


75 ± 31 ± 26 


WA92 (98) 


m 






f- 7+1.3+0.6 
J- ' -1.1-0.5 


E789 (95) 


m 


< xp < 0.1, < 2 




5.7±1.5±1.3 


E771 (99) 


m 


-0.35 < x'l^^' < 0.15 




43l^^ ± 7 


HERA-B (06) 


m 


13.3 ±2.9 


14.9 ±2.2 ±2.4 



Table 9: BB pair production cross-sections measured in fixed-target experiments, in 
the probed phase space (A) and extrapolated to full phase space (4-7r). NAlO quoted 
a second value, 10 ± 5 nb, obtained with different assumptions on the B meson 
kinematical distributions. 
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Figure 11: Fixed-target BB cross-section measurements, including (open square) the 
cross-section NAIO derived by using alternative B hadron kinematical distributions. 



a different decay channel for each of its 18 B mesons. . . Also WA92 explored several 
different event topologies, taking the branching ratios from the PDG 1996 tables |75j . 

In Table El and Fig. ^2 we summarise the BB pair production cross-sections mea- 
sured in fixed-target experiments, as taken from the original publications. Like in 
the charm sector, all experiments assumed a linear scaling of the beauty production 
cross-section with the mass number of the target nucleus to extract the elemen- 
tary hadron-nucleon production cross-section, irrespective of the phase space window 
where the measurement was made. 
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We will now describe the models used to evaluate detector acceptances and effi- 
ciencies, and to extrapolate the measured cross-sections to the full phase space values 
compiled in Table IHl and Fig. ^2 We also mention some selection cuts applied on 
the kinematics of the detected decay products, often not easy to relate with the B 
mesons' kinematical distributions. 

• NAIO (88) The NAlO Collaboration measured the BB cross-section by trig- 
gering on dimuons and looking for a third muon in the offiine analysis. This 
decay pattern is obtained when both B mesons decay through the semi-muonic 
channel, B — > DfiX, and one of the D mesons also decays into a muon. An 
offiine cut of pi^ > 1.4 GeV/c was required for two out of the three muons. At 
the times of NAIO (the article was submitted in 1987) the theory of beauty 
production was not yet developed in detail and a simple model was used with 
the following characteristics: the BB pair is assumed to be the only decay prod- 
uct of an intermediate state with the mass of the T'". This state is produced 
at mid-rapidity with a gaussian distribution in mass, with a width depending 
on the centre-of-mass energy [TH]. The kinematical distributions were taken to 
be d^cr/dxpdpT oc (1 — |xf|)^Pt exp(— 2pt)- The general features concerning 
the B and D meson decays were taken from Ref. [77j. In order to compute 
the BB cross-section, beauty branching ratios were taken from the PDG 1986 
tables Uni; the value BR(D /i + A) = 11.0% was deduced from Ref. [HI, 
assuming that D mesons are 75 % neutral and 25 % charged. With these assump- 
tions, a BB cross-section of 14^*^5 nb/nucleon at Ei^h = 286 GeV was obtained. 
Alternatively, when NAIO used the x-p and distributions of Ref. [7H], a BB 
cross-section of 10 ± 5 nb/nucleon is derived. In both cases only the statistical 
error is given. 

• WA78 (89) The WA78 Collaboration explored the same decay modes as NAIO, 
71^ + U ^ BBA with B — > D/iA and D — /iA, triggering on two or more 
muons. 68 BB candidates (including an expected background of 5.2 events) 
were found in the "like-sign" dimuon events passing the following requirements: 

Evis = Ecal + EE^< 300 GeV, |>T,tot = EPT > 2.7 GeV/c, Eiept = ^beam - 

Ecal > 100 GeV, where Ecai is the total energy deposited in the calorimeter. 
Another 11 candidates (including 1.1 estimated background events) were found 
in the "three muon sample", with E^is < 270 GeV and px.tot > 3 GeV/c. Two 
theoretical models were used to extract the BB production cross- sect ion. The 
ffist one was a "LO QCD model", based on the calculations outlined in Ref. j79j . 
with d^cr/dxFidxF2dpT = d^cr/dxFidxF2 exp(— p^/c), where xfi and xf2 are 
the Feynman x of the b and b quarks and pt their transverse momentum, 
with c ~ 6.9 (GeV/c)^. From the 2-dimensional distribution d^(T/da:Fida;F2 
given by the model, the single b quark distribution could be represented by 
dcr/dxF oc exp(— (xf — 0.09)^/^4^), with ^ 0.3. In order to evaluate the 
sensitivity of the results to the mean value of xp, its value was changed from 
(xf) = 0.09 to 0.05, increasing the cross-section value by 25 %. In order to 
compare the WA78 results to those of NAIO, also the "NAIO model", described 
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above, was used to extract the BB cross-section, yielding a value higher by 
a factor of ~3.8, due to the very different xp distributions assumed in the 
two models. The B meson decays were simulated with the Lund Monte Carlo 
program. If WA78 uses the "NAlO model" , the two measurements are in good 
agreement, but when the QCD-based production model is used, which seems to 
provide a better description of the observed Eyept distributions, the agreement 
suffers significantly. The B°'s in the observed event sample were an unknown 
mixture of Bj^ and Bg . Apart from using different production models, also three 
different values for the mixing parameter were explored, xb = 0, 0.1, 0.2. In 
the PDG 2004 edition the value for the mixing parameter for an unknown 
mixture of B meson species is Xb = 0.1257 ± 0.0042. Therefore, out of the 8 
different values given in the publication for the BB production cross-section, 
we retained the one corresponding to Xb — ^-^ LO QCD model, 3.6 ± 

0.4 ± 1.1 nb/nucleon. The systematic error, which is higher than the statistical 
one, includes uncertainties on the acceptance and track reconstruction efficiency 
calculations, on the absolute normalisation and on the branching ratios. 

• E653 (93) The E653 Collaboration searched for beauty candidates by looking 
at topological events with muonic secondary vertices and an associated high 
Pt muon (px > 1.5 GeV/c). The and xp distributions of the reconstructed 
18 beauty mesons were compared to LO [ZHl EDI and NLO HI El ESI QCD 
calculations, and to Pythia jHlj. Within statistical errors, all models gave a 
satisfactory description of the experimental data and their differences resulted 
in a 10 % systematic uncertainty in the shape of the differential cross-section. 
Other sources of systematic uncertainties include the luminosity calculation 
(5%), the semi-muonic branching ratio (8%), the detection efficiency model 
(7%), and others (10%). 

• E672/E706 (95) The E672/E706 Collaboration studied the B ^ J/t/'X de- 
cay. The evaluation of acceptances and detection efficiencies of beauty hadron 
pairs was done using the MNR NLO calculations [S31EH]! with hadron momenta 
equal to that of the parent quark, i.e. without smearing of the momenta due 
to fragmentation. In their phase space window, xp > 0, they measured a BB 
cross-section of 47 ± 19 ± 14 nb/nucleon, based on 8 ± 3.3 signal events and a 
background of 2 ± 1. The systematic error includes uncertainties in the normal- 
isation (13%), branching ratios (13%), b quark production, hadronisation and 
decay properties (16%), and reconstruction efficiencies (18%). The extrapola- 
tion to full phase space was done using the MNR NLO calculations jHS] but the 
systematic error also reflects the use of alternative xp distributions |38^ ITH] in 
the extrapolation. 

• WA92 (98) The WA92 Collaboration classified the beauty candidates in three 
samples, according to the event topologies. The first class, composed of semi- 
muonic decays, had 12 candidates. The second class collected 12 multi-vertex 
events, given the large decay multiplicity of beauty events. In the third class, 
the decay channel B — ^ DX was investigated, giving 5 candidates. To evaluate 
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the acceptances, Pythia [Sni and Jetset jSTj were used with the EHLQ [HE] 
parton distributions for the nucleons in the target and the Owens jHHI PDFs for 
the incoming pions. The b quark fragmentation was simulated by Jetset, using 
the Lund string fragmentation model fH]- The PDG 1996 [75] branching ratios 
were used. The BB cross-section was then obtained from the three event samples 
individually, yielding compatible results: = 6.2t.l'j'tQ% (Tmvtx = 5.2^^;3lQ;g 
and (Td = B.Olg'g^o s nb/nucleon. Combining these values gave an overall full 
phase space beauty cross-section of 5.71};^1q;5 nb/nucleon. The systematic error 
includes a 6 % uncertainty from the luminosity calculation and uncertainties due 
to the acceptance and efficiency calculations, which were 9%, 12% and 17% 
for the first, second and third event samples, respectively. The uncertainties in 
the branching ratios were not taken into account. 

• E789 (95) For the evaluation of acceptances and efficiencies, E789 took the 
xp and pt distributions of the b quarks from the MNR NLO model jHOl IHHj : 
the intrinsic transverse momentum was simulated with a Gaussian of sigma 
0.5 GeV/c. The b quark fragmentation was modeled with the Peterson frag- 
mentation function with e = 0.006 ± 0.002 jUJ; and the average life- 
time of B hadrons was set to 1.537 ± 0.021 ps, as given in the PDG 1994 
tables [70], which also provided the values of the branching ratios. The J/^/^'s 
coming from B — * J/tpX decays were simulated with the momentum and de- 
cay angle distributions measured by the CLEO Collaboration. The measure- 
ment was carried out in the < Xp'^ < 0.1, p^^ < 2 GeV/c kinematical 
region, representing 15 % of the full phase space (value quoted in the original 
publication and reproduced by Pythia). Within this phase space domain, a 
cross-section of (d^cx/dxp dp^) = 81 ± 21 ± 15 pb/(GeV/c)^/nucleon was mea- 
sured. The extrapolation to full phase space yielded a total cross-section of 
5.7±1.5ibl.3 nb/nucleon. The systematic error is dominated by the evaluation 
of the luminosity (11%), efficiencies (10%), b quark production, hadronisation 
and decay models (8%), fit of the mass spectrum (5%), and J/tp branching 
ratio (4%). _ 

• E771 (99) The E771 Collaboration measured the BB production cross-section 
mainly from the semi-muonic decay of both B mesons, BB — > /i"*"/!". Events 
in which one beauty hadron decayed through a J/ip, B J/ipX fi~^fi~ X, 
also contributed to the evaluation of the beauty cross-section. The muons were 
detected in the range —0.25 < Xp < +0.5, and in order to pass the muon filter 
they needed to have p^^ > 6 GeV/c (or p^ > 10 GeV/c, in the angular region near 
the beam). We estimated, using Pythia 6.208, that this phase space window 
covers 88 % of full phase space. E771 evaluated acceptances and efficiencies 
using Pythia 5.6 and Jetset 7.3 [021 • -^^^ semi-muonic decays, BB — > fi'^fi', 
the default Pythia branching ratios were taken. The main sources of systematic 
uncertainties concern the production model (10%), the luminosity evaluation 
(5%), the efficiencies (9%), and the semi-muonic branching ratio (7%). 

• HERA-B (06) To minimise the systematic uncertainties due to luminosity and 



32 



efficiencies, HERA-B measured the BB production cross-section relative to the 
prompt J/V' production cross-section. The J/V' mesons due to decays of beauty 
hadrons were distinguished from those directly produced at the target thanks 
to the good vertex resolution, as done by other experiments. The b-content of 
the selected sample was con&med by a lifetime measurement, based on the ob- 
served decay length. The beauty events were generated with Pythia 5.7 [^21 and 
weighted with probability distributions obtained from NLO calculations [5Hll81j . 
using NNLL MRST PDFs g Ul, a b qu ark mass = 4.75 GeV/c^, and a QCD 
renormalisation scale /i = ^/m'^ + p^. The colliding partons were assigned an 
intrinsic transverse momentum distributed according to a Gaussian with sigma 
0.5 GeV/c. The fragmentation was described by the Peterson function, with 
e = 0.006. The subsequent B hadron production and decay were controlled 
by Pythia's default settings. The "prompt" J/ip events were generated tak- 
ing into account the differential cross-sections, da/dp"^ and da/dxp, measured 
by E771 in p-Si collisions, to properly consider nuclear effects (such as the 
"Cronin effect"). The reference J/ip production cross-section, at the energy of 
HERA-B, was derived from a global analysis of all published J/ijj cross-section 
measurements, including the value measured by HERA-B, with the help of 
non-relativistic QCD calculations. A total beauty production cross-section of 
14.9 ± 2.2 ± 2.4 nb/nucleon was obtained by combining the analyses in the 
muonic and electronic channels of all collected data, and by extrapolating to 
full phase space the measurement made in the —0.35 < xp < +0.15 window 
(covering 90.6±0.5 % of Air). The overall systematic uncertainty is 14 %, mainly 
due to the B — > J/ipX branching ratio (8.6%), the trigger and reconstruction 
efficiencies (5 %), the b production and decay model (5 %), and the prompt J/ip 
production model (3.1%). 

At collider energies, two experiments, UAl and CDF, measured beauty produc- 
tion cross-sections in pp collisions. Before we mention specific numerical values, we 
should clarify that the experiments are not consistent in the way they quote their 
measurements, referring to the production cross-sections of beauty mesons, beauty 
hadrons or beauty quarks, sometimes meaning the single flavour cross-section (only 
beauty, not anti-beauty), etc. Each paper accurately describes what was measured, 
but it is not trivial to get a consistent picture of all the measurements. 

The UAl Collaboration measured the beauty production cross-section at a/s = 
630 GeV, by combining four independent analyses jHj . For each of these analyses 
the systematic error was evaluated including uncertainties from the luminosity evalu- 
ation, acceptance and efficiency calculations, background estimation, and additional 
analysis-specific uncertainties. In order to relate the measured cross-sections to the 
production cross-sections of beauty hadrons, the Monte Carlo model ISAJET jUlI was 
used, as described in detail in Ref. [HH]. The Peterson fragmentation function was 
used, with e = 0.02 (softer than the standard 0.006 value). Systematic errors due 
to this deconvolution include uncertainties on the fragmentation function (6%), on 
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the branching ratios and on the assumed shape of the b quark's transverse momen- 
tum (20%). The inclusive b quark and B hadron cross-sections for \y\ < 1.5 and 
Pt > PT,min, obtained in this way for each individual analysis, are summarised in 
Table 2 of Ref. jlj. The combined cross-section value, extracted for p\ > Q GeV/c 
and \y\ < 1.5, was then extrapolated to full phase space using a NLO QCD calcula- 
tion jnOlEI], with factorisation scale /i = -\/ (pT.min)^ + and mb = 4.75 GeV/c^. 
The beauty cross-section extrapolated to PT,min = was dominated by the "low-mass 
dimuon" and "J/V'" analyses, which used pf^ > 3, pt^ > 2, and pl^^ > 5 GeV/c. It 
gave a b quark (beauty only, not anti-beauty) production cross-section in the central 
rapidity range, \y\ < 1.5, of cr(b) = 12.8 ± 4.7(exp) ± 6(th) /ib. Using the rapidity 
dependence predicted by the model, UAl derived a total beauty production cross- 
section of (T(bb) = 19.3±7(exp) ±9(th) /ib. The first error is due to the normalisation 
of the theoretical QCD shape to the data; the second error is due to the extrapolation 
to full rapidity, from uncertainties on the shape of the QCD curve. 

In Run I, CDF obtained the beauty production cross-section at y/s =1.8 TeV 
by reconstructing the exclusive decay — > J/ipK — > /i^/z" K, selecting kaons of 
Pt > 1.25 GeV/c. For the evaluation of the acceptances, a Monte Carlo simulation 
based on a NLO QCD calculations [SIIIH2] was performed, using the MRST PDFs jHS] 
and mb = 4.75 GeV/c^. The renormalisation and fragmentation scales were both set 
to -\/m^ -|- The hadronisation into B mesons was modeled with the Peterson 
fragmentation function, with e = 0.006. The B J/ipK decays were simulated 
using a modified version of the CLEO Monte Carlo program [HZj. The B"*" meson 
production cross-section was measured in the phase space domain p^ > 6.0 GeV/c 
and \y\ < 1.0, and was not extrapolated to full phase space. The value measured 
was (Jb+{Pt > 6.0, \y\ < 1.0) = 3.6 ± 0.4(stat © systp^) ± 0.4(syst) /xb, the first error 
being the quadratic sum of the statistical and px-dependent systematic errors, while 
the second is the same for all bins. The main contribution to the systematic 
error is the uncertainty on the (combined) branching ratio (10.2%). Other sources 
of systematic uncertainties include the luminosity evaluation, the influence of in- 
flight kaon decays, and pT-dependent trigger efficiencies. From the theoretical side, 
Pt dependent uncertainties were estimated due to the QCD renormalisation scale 
and to the Peterson parameter, e. Meanwhile, the uncertainty on the combined 
branching ratio decreased from 10.2 % to 4.3 %, so that in this case an update would 
decrease the overall systematic error (see Table |H1). Note that the published value is 
the B"^ production cross-section, obtained dividing by 2 the sum of the B"*" and B~ 
measurements. The production cross-section of beauty hadrons, of a single flavour 
(only beauty or only anti-beauty) would be a factor 2.5 higher, according to the 
PDG 2004 tables [021, where the B"*" appears as being 39.7% of all beauty hadrons 
produced, the rest being B° (39.7%), B° (10.7%), and b-baryons (9.9%). 

In Run II, at ^/s = 1.96 TeV, CDF searched for beauty decay topologies in 
299 800 events with a reconstructed J/tp, leading to 38078 beauty hadrons. For the 
acceptance evaluations, the simulated distributions (77, px, ;2-vertex, offset, . . . ) were 
compared to the reconstructed data and the input distributions tuned until there 
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was perfect agreement between data and MC. For cases where the quahty of the data 
is good enough, such an iterative procedure is the best possible method to evaluate 
acceptances, since it does not rely on theoretical inputs. Furthermore, since CDF did 
not extrapolate the measured cross-section to full phase space, no theoretical model 
was used at all. The systematic errors are, hence, entirely of experimental origin. 
Very detailed descriptions of the ^/ip and beauty analyses are given in Ref. jSS]- The 
inclusive single B hadron cross- section is a(pp ^ BX, \y\ < 0.6) = 17.6 ± 0.4^1^ /ib. 
To be precise, this is the cross-section for one B hadron to be within the \y\ < 0.6 
window, irrespective of where the second B is. We emphasise that this is the cross- 
section for the production of any beauty hadron (mesons and baryons; charged or 
neutral), but only beauty or only anti-beauty, not both flavours. Hence, this is the 
same numerical value as the BB pair cross-section, (t(BB). It is also the same as the 
b quark pair cross- sect ion, cr(bb), since bottomonium production is negligible. 

CDF observed that the Run I measurement of the B''" production cross- 
section, for > 6.0 GeV/c, was higher than the NLO calculations |HI1 IH3, by a 
factor 2.9 ± 0.2(stat © systp^) ± 0.4(syst). However, Cacciari and Nason jHEl pointed 
out that the discrepancy is not really due to the perturbative calculations of the cross- 
section but rather to the use of the Peterson fragmentation function with e = 0.006. 
Using smaller values of e in the framework of a NLL (next-to-leading-log) calculation 
reduced the discrepancy to 1.7 ± 0.5 (experiment) ± 0.5(theory) (including a better 
treatment of the theoretical error bars). For more information on this topic, see 
Refs. [HHt inn] • It is worth mentioning that the Run II and Run I measurements agree 
with each other, when compared for > 6.0 GeV/c. 

To clarify the source of the disagreement between the CDF data and the NLO 
QCD calculation, the Tevatron ran for nine days at ^/s = 630 GeV |100j . so that 
CDF could measure the ratio between the beauty quark production cross-sections at 
630 and 1800 GeV. The UAl measurement, at 630 GeV, had not shown a significant 
departure from expectations, but was affected by a larger uncertainty. It is obvious 
that, both experimentally and theoretically, it is much more accurate to obtain a 
ratio of cross-sections at two different energies than each of the absolute values, since 
many uncertainties cancel out. Given the very short beam time available, a highly 
efficient beauty finding algorithm had to be developed, using single moderate-px 
muon triggers, combining the muon with a high-px track and selecting the high- 
mass muon-track combinations with positive lifetime. The two data sets (630 and 
1800 GeV) were as similar as possible, with all data collected between December 1995 
and February 1996, and with identical online and offline event selections. The cross- 
section ratio was reported for px of the beauty quark above 10.75 GeV/c, because 
Monte Carlo simulations showed that 90 % of the reconstructed and identified beauty 
events have a px above this threshold. The measured ratio, between the beauty 
quark production cross- sect ions for px above 10.75 GeV/c at 630 and 1800 GeV, is 
0.171 ± 0.024(stat) ± 0.012(syst), in very good agreement with the theoretical NLO 
QCD prediction jHU 1101, shown in Figs. 3 and 4 of Ref. |100j . The conclusion 
is that the disagreement between the CDF data and the NLO QCD calculation is 
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essentially the same at 630 GeV and at 1800 GeV. It can also be seen that the CDF 
beauty cross-section at 630 GeV, derived by combining the ratio with the absolute 
value measured at 1800 GeV, and extrapolated from \y\ < 0.6 to \y\ < 1.5, is about 
a factor 2 higher than the value of UAl, for roughly the same pt range, as shown in 
Fig. 5 of Ref. fTUn] . 

6 Pythia calculations versus data 

After having collected and revised the open charm and open beauty production cross- 
sections measured by many experiments, we will now compare them with theoretical 
calculations made using the Pythia Monte Carlo event generator |llj . 

Our first motivation is to describe the energy dependence of the measured pro- 
duction cross-sections, so as to calculate the values expected for energies relevant to 
recent or on-going experiments at the CERN SPS (NA60) and at RHIC (PHENIX 
and STAR). However, the yield of D or B mesons to be expected in a given detector 
depends not only on the production cross-section but also on the kinematical distri- 
butions of the single mesons, on the pair correlations, etc. Therefore, we will also see 
in this section how the results of Pythia compare to the available differential data, 
and how sensitive they are to changes in certain settings. 

Pythia is an easily accessible tool which has been and continues to be extensively 
used by many experiments interested in heavy flavour production (to evaluate accep- 
tances, efficiencies and other important elements of the data analyses procedures). It 
may very well be, however, that better calculations of heavy flavour production, in 
one or another aspect, are provided by other computer codes, such as Herwig |101j . 
IS A JET inU, MC@NLO [1021, or MNR NLO |HS1EI]. Surely, a work similar to the 
one we present in this section could and should be made using those other theoretical 
models. 

In some of our calculations, we varied the set of PDFs, the mass of the charm quark 
and the deflnition of the scale, as will be explained below. The intrinsic transverse 
momentum of the colliding partons, k^, was generated according to a Gaussian distri- 
bution with a width determined by the parameter PARP(91): exp(— /c^/PARP(91)^). 
This expression is often written as exp 

(-fc2/(A;2))^ deflning PARP(91) as ^/i^). We 
used the Lund string fragmentation scheme ^H] ? modifled for heavy flavoured quarks 
by Bowler ^7j, which provides a good description of the high-precision data obtained 
by the SLD, ALEPH, DELPHI and OPAL Collaborations [TI7n| . 

We generated charm and beauty events setting MSEL = 4 and 5, respectively, 
to ensure that every generated event has a cc (or bb) pair. This is particularly 
important when simulating collisions at relatively low energies, where heavy flavour 
production is a rare process. Besides the LO diagrams (quark-antiquark annihilation 
and gluon fusion, see Section |2I), the calculations also include initial and flnal state 
radiation. However, they do not include "flavour excitation" and "gluon splitting" 
diagrams. These can also be generated by Pythia, but only in the context of a 
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full minimum bias generation (MSEL = 1), a very inefficient operation mode at 
low collision energies (see Refs. jl(J4l I1U5] for details). The absence of higher-order 
diagrams is (partially) compensated with scaling up the calculations by empirical 
K-factors, under the hypothesis that the kinematical distributions do not change 
significantly between leading and higher-order diagrams. 

6.1 Absolute cc production cross-sections 

As a ffist step, we used Pythia's default values for the mass of the c quark {rric = 
1.5 GeV/c^) and for the definition of the scale, while varying the set of parton dis- 
tribution functions. For the proton PDFs we used the CTEQ6L (2002) |], MRST LO 
(2001) [lU] and GRV LO (1998) sets. Since NLO PDFs have also been used 

with Pythia in the past, despite the inconsistency resulting from the fact that Pythia 
does not include higher-order diagrams, we have also used the CTEQ6M (2002) and 
MRST c-g (2001) sets for the extraction of the K-factors. To describe the parton 
densities in pions we used the following sets: GRV-P LO (1992), SMRS-P2 (1992), 
GRV-P HO (1992) and ABFKW-P Set 1 (1989). As indicated by the year, the avail- 
able pion PDFs are considerably older than the proton PDFs. Furthermore, there 
is only one set of LO pion PDFs. For the simulations with it beams we have used 
MRST LO (2001) to describe the target nucleons, since they were fit with a Aqcd 
value, 220 MeV, similar to the values used in the vr PDFs, Aqcd = 190-231 MeV. 

To easily compare the results obtained, we parameterised the y/s dependence of 
the cc production cross-section with the expression 



inspired by a formula commonly used to describe the energy dependence of J/ip 
production |l()7j . The parameter ps was fixed to 0.35 in pp collisions and to 0.3 in 
vr^p colhsions, while the other parameters were adjusted to each specific set of PDFs. 
The calculated curves for neutral and charged D meson production in pp and 7r~p 
collisions are presented in Fig. where they are compared to the corresponding 
experimental measurements. Tables ID and El normalised to xp > 0. 

In several cases, the same measurement is represented twice, before (open marker) 
and after (solid marker) the corrections explained in Section |3] Note that the mea- 
surements of E743, E653 and E789 were all performed at -Eiab = 800 GeV; to improve 
their visibility, they were slightly displaced in Fig. (without affecting any calcula- 
tions). The same is valid for the NA16 and NA27 7r~p measurements, both performed 
at 360 GeV. Each curve appears twice in the figures, as directly calculated by Pythia 
and after being scaled up by a K-factor fitted to the data. The solid lines corre- 
spond to the LO PDF sets. The E789 p-Au measurement of the neutral D meson 
cross-section was not considered in the fitting of the K-factors, since it would have 
significantly deteriorated the quality of the fits, maybe because this measurement was 
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Figure 12: Energy dependence of the production cross-section of neutral (left) and 
charged (right) D mesons, in proton (top) and pion (bottom) induced collisions. The 
upper curves were scaled by K-factors fitted to the data (excluding the E789 point). 

made in an exceptionally small phase space window (0 < xf < 0.08, pt < 1-1 GeV/c). 
We remind that the HERA-B measurement is still preliminary. 

Using the appropriate K-factor for each curve, we can see from Fig. that all of 
them describe the ^/s dependence of the data points within the experimental error 
bars. The reduced values of the fits are around 1 or lower for the proton data and 
around 2 for the pion data. 

Table ^1 shows the K-factors which best describe the four sets of data. To re- 
produce the measured data points, the K-factor required by the charged D mesons is 
significantly higher than the one needed by the neutral mesons. For instance, in pp 
collisions the CTEQ6L and MRST LO sets require K-factors around 5 and 2.5 for the 
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Table 10: K-factors obtained from fitting the experimental data points with Pythia's 
default settings, for several sets of PDFs. The K-factors of the proton and pion data 
have a relative uncertainty of 9 and 5 %, respectively. The MRST LO (2001) PDFs 
were used for the target nucleons in the vr" induced collisions. 



charged and neutral D mesons, respectively. The significant difference between the 
charged and neutral D meson K-factors, required by the experimental data, indicates 
a problem in the probabilities used by Pythia when fragmenting the charm quarks 
into D hadrons: in ~ 63 % of the cases, Pythia fragments the charm quarks into D° 
or DO mesons, while the probability of producing D"*" or D~ mesons is only 20 % (the 
remaining charm quarks hadronise essentially into Dg or A^). This is due to the larger 
feed-down contribution from D* decays to neutral D mesons. It is easy to derive the 
ratio between the charged and neutral D meson yields as given by Pythia; it follows 
directly from the assumption that the c quark has equal probabilities of fragmenting 
into the neutral and charged D (singlet) or D* (triplet) states, and from the numerical 
values of the branching ratios of the D* to D decays [S^ : 

a(D+) _ ajDj^^J + a(D-+) ■ B{D*^ ^ D+) + a(D-°) ■ BjD*^ ^ D+) _ 

a(DO) a(DO_J + a(D*+) ■ B{D*+ D^) + a(D*o) . B{D*o D^) 

0.25 + 0.75 -0.323 + 0.75 -0.0 ^ ^ 

= ! . = 33 

0.25 + 0.75 -0.677 + 0.75- 1.0 

Figure El shows that this "expected value" is considerably lower than the ratios 
given by the hadro-production experiments which measured both charged and neutral 
D mesons. The average values of the measured ratios are 0.61+0.08 in p-A, 0.37+0.03 
in 7r~p collisions, and 0.41 ± 0.03 if we merge both data sets. This discrepancy is 
due, at least in part, to the fact that (naive) spin counting is not applicable to 
charm production: it would imply that the charm quark fragments into a vector 
state with a probability Py = 3/(3 + 1) = 0.75, while all available measurements 
(in hadro-production, photo-production, at LEP, etc) give a combined value of Py = 
0.59 ± 0.01 |108j . as illustrated in Fig. El This is presumably because the mass 
difference between the D and D* mesons cannot be neglected, making the lighter D 
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Figure 13: Ratio between charged and neutral D meson cross-sections, for pp (left) 
and 7r~p (right) collisions, compared to the values 0.33 and 0.42, expected when using 
Pv = 0.75 and 0.6, respectively. 
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Figure 14: Fraction of directly produced vector D* states to the total, vector and 
pseudo-scalar states, averaged over the measurements made for each of the indicated 
production systems. The central band represents the world average 0.596 ± 0.012, 
quite different from the spin-counting expectation, 0.75. Figure taken from Ref. |108j . 

mesons more probable to be produced. In Pythia, the value of Py is given by the 
parameter PARJ(13), which is 0.75 by default. Setting it to 0.6 provides a much 
better agreement with all the presently available measurements, as just mentioned, 
and increases the ratio between the charged and neutral D meson production cross- 
sections from 0.33 to 0.42. 

We will now assume that the ratio between charged and neutral D meson pro- 
duction cross-sections is the same in proton and pion induced collisions, as suggested 
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by the Py universality seen in Fig. (and consistent with the left panel of Fig. IT^ 
given the large error bars of the proton data). This means that we can use a common 
K-factor for charged and neutral D mesons, as long as we set Py = 0.6. Table ^2 
gives the K-factors extracted in this way, for several sets of PDFs. 



PDF set 


K-factor 


CTEQ6L (2002) 


3.0 


MRST LO (2001) 


3.8 


GRV LO (1998) 


4.6 


CTEQ6M (2002) 


2.6 


MRST c-g (2001) 


2.8 



Table 11: K-factors obtained from simultaneously fitting the charged and neutral D 
meson p-A data with Pythia calculations, setting PARJ(13) = 0.6. These values have 
an uncertainty, from the fit, of ~ 8 %. 

It is unfortunate that, until now, Pythia uses the same parameter, PARJ(13), to 
also define the Py value for beauty production, where the value 0.75 seems to be 
appropriate. Indeed, measurements made by ALEPH, DELPHI, OPAL and L3 give 
an average value of Py = 0.75 ± 0.04 |108j . It would be more useful to have different 
parameters for charm and for beauty, so that they could be given different values (as 
done for the u/d and s quarks, through the PARJ(ll) and PARJ(12) parameters). 
This is not so crucial if the charm and beauty events are independently generated 
(with MSEL = 4 or 5, respectively) but becomes mandatory if we generate both 
heavy flavours simultaneously, as frequently done. 
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Figure 15: xp distributions for charged and neutral D mesons, as measured by WA92 
in pion induced collisions. The inset shows their ratio. 

FigureElshows that the ratio between charged and neutral D meson cross-sections 
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does not seem to depend on xp, according to the WA92 measurements [IH]- This ratio 
also seems to be independent of the colhsion energy, when comparing measurements 
made by many experiments (see Fig. [T3j). 

We will now concentrate on the total cc production cross-sections for pp collisions, 
in our evaluation of the influence of certain parameters on Pythia's results. 




Figure 16: Total cc production cross-sections for fixed-target energies (left) and up 
to ^/s = 200 GeV (right). Open symbols indicate the pp measurements. 







PDF set 


K-factor 


xVndf 


^lab = 158 


^lab = 400 


v/i= 200 


CTEQ6L (2002) 


3.0 


1.4 


3.6 


17.5 


803 


MRST LO (2001) 


3.8 


0.7 


4.8 


18.4 


439 


GRV LO (1998) 


4.5 


0.7 


5.2 


18.3 


563 


CTEQ6M (2002) 


2.5 


0.7 


4.6 


18.4 


427 


MRST c-g (2001) 


2.7 


1.0 


4.0 


18.0 


524 



Table 12: K-factors which provide the best description of the cc data in pp and p-A 
collisions, for each PDF set. The values have a relative uncertainty of around 7%. 
The last three columns give the elementary pp cross- sections calculated by Pythia 
with these K-factors, for three different energies, given in GeV. 



FigurelTI)! compares the measured cc production cross-sections with Pythia's curves, 
obtained with five different PDF sets. The curves are labelled in the order, from top 
to bottom, in which they occur at ^/s = 15 and 200 GeV, in the left and right panels, 
respectively. The extracted K-factors for each PDF set are summarised in Table IT^ 
together with the corresponding pp cross-sections for typical SPS and RHIC energies. 
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At a/s = 200 GeV, the calculated pp cross-sections vary between 440 and 800 /ib, if 
we only consider LO PDFs. 




Figure 17: LO cc production cross-sections for rric = 1.3, 1.5 and 1.7 GeV/c^, scaled 
up by the appropriate K-factor. The inset shows the corresponding ratios. See the 
text for details. 

In a second step we studied the influence of varying the mass of the c quark by 
~ 15 % with respect to Pythia's default. We used the CTEQ6L PDFs for these calcu- 
lations, but other proton PDFs give similar results. Figure IT7I shows the calculated cc 
production cross-section in pp collisions, using rric = 1.3, 1.5 and 1.7 GeV/c^, scaled 
up by the appropriate K-factor. While smaller rric values lead to higher calculated 
cross-sections, in particular at energies close to threshold, after the curves are nor- 
malised using the available fixed target data, and given the somewhat different shapes 
of the calculated curves, it turns out that at = 200 GeV the estimated cc cross- 
section is 35 % higher with rric = 1-7 GeV/c^ and 30 % lower with rric = 1-3 GeV/c^, 
with respect to the default value. The results are summarised in Table ITSl 





(Tec [/ib] 


rric [GeV/c2] 


K-factor 


xVndf 


^lab = 158 


^lab = 400 


v/i = 200 


1.3 


1.2 


0.9 


4.5 


18.1 


569 


1.5 


3.0 


1.4 


3.7 


17.5 


811 


1.7 


6.6 


2.0 


3.2 


16.9 


1100 



Table 13: Same as previous table, when varying the mass of the c quark, rric- The 
CTEQ6L PDFs are used. See the text for details. 



Different definitions of the squared energy- momentum transfer, Q^, can be used. 
To evaluate the influence of this setting on our results, we replaced Pythia's default. 
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equivalent to = in the processes we are studying, by = s, the choice of 
Refs. |1(J9U11U] . Figure ITHl shows the effect of using these two different definitions 
on the cc cross-section, keeping rric = 1.5 GeV/c^ and using the CTEQ6L PDFs. 




Figure 18: cc production cross-sections, for two different definitions, scaled up by 
the appropriate K-factors. The inset shows the ratio of the curves, taking Pythia's 
default as the reference. 

We see that using s as the definition leads to significantly lower cross-sections 
with respect to the values obtained when using Pythia's default setting. The dif- 
ference is energy dependent: at low energies the cross-sections obtained with the s 
definition are around 3 times lower, while at ^/s = 200 GeV the difference reduces to 
a factor of 2. Once the curves are scaled up to describe the data, the steeper rise with 
^/s of the = s curve leads to 60 % higher cross-sections at ^/s = 200 GeV, with 
respect to the values obtained with the default setting. The results are summarised 
in Table El Calculations with other PDF sets give comparable results. 







definition 


K-factor 


xVndf 


^lab = 158 


^lab = 400 


v^ = 200 


default 


3.0 


1.4 


3.8 


17.5 


808 


s 


8.8 


1.9 


3.4 


16.9 


1327 



Table 14: Same as previous tables, when varying the definition, keeping rric = 
1.5 GeV/c^ and using the CTEQ6L PDFs. Note the high K-factor required by the 
Q"^ = s choice. 



These calculations show that the cc production cross-section at ^/s = 200 GeV, 
as derived from Pythia's calculations normalised by the existing SPS, FNAL and 
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HERA-B measurements, can vary by ±30 % due to the use of different sets of PDFs 
and by around ±30% if tlie c quark mass is changed by ±15%. Furthermore, using 
= s, as done by some experiments, leads to a 60 % higher cc cross-section at 
^/s = 200 GeV. From Table ^1 where we used the default c quark mass and 
definition, and only considering the values corresponding to the LO PDF sets, we 
obtain 600 ± 30 % /xb as our best estimate for the cc cross-section at ^/s = 200 GeV. 
The uncertainties of the calculations for SPS energies are smaller. At i?iab = 158 
and 400 GeV, we expect total cc cross-sections of 4.5 ± 20 % /xb and 18 ± 5 % fib, 
respectively. 

6.2 Single D meson kinematical distributions 

In this section we compare measured D meson kinematical distributions to calcula- 
tions done with Pythia. These distributions are sensitive to non-perturbative effects, 
such as the intrinsic transverse momentum which the partons carry before they collide 
and the fragmentation of the heavy quarks into hadrons. 

Table ITHl gives an overview of the experiments which measured pt and x-p distribu- 
tions of D mesons produced in p-A and vr-A collisions. Some experiments separately 
give the distributions of "leading" and "non-leading" particles. A D meson is called 
"leading" if its c or c quark combines with one of the non-interacting beam or tar- 
get valence quarks. These quarks have in general a much higher momentum than a 
light quark from the sea, so that the formed D meson will have a visible (forward 
or backward) boost. In pp and p-A reactions the D~ and are the leading parti- 
cles. In 7T~ induced reactions, for xp > 0, the D~ and D° are leading. This purely 
non-perturbative (hadronisation) effect is properly described by the Lund string frag- 
mentation model, as shown in Refs. |104t I105j . See end of next section for further 
information. 

We will now focus on the and xp distributions of single D mesons. The E791 
Collaboration measured the and xp distributions of neutral D mesons, produced 
in 500 GeV 7r~-C collisions jlU]- This measurement represents by far the largest D 
meson data sample at fixed-target energies, with ~ 90 000 fully reconstructed neutral 
D mesons. Figures 5 and 6 of Ref. jlOj show comparisons between the measured 
distributions and curves calculated with Pythia 5.7 and Jetset 7.4 jHSj- With these 
code versions, of 1994, the calculations were unable to reproduce the measured dis- 
tributions. On the contrary, the most recent version of Pythia, 6.326, describes quite 
well ^ the a.nd xp distributions of E791, with PARP(91) between 1 and 1.5 GeV/c, 
as can be seen in Fig. ^| 

Note that from version 6.314 (Oct. 2004) onwards the default value of PARP(91) 
changed from 1.0 to 2.0 GeV/c. The E791 calculations, made with version 5.7, used 

^When we made calculations with version 6.325 we noticed that the resulting distributions 
were abnormally hard. T. Sjostrand, main author of Pythia, immediately identified the problem, 
introduced in version 6.319, and sent us the solution, later on implemented in Pythia 6.326. 
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Experiment 


Beam i?iab 


Target 


particle 




events 








[GeV] 










NA27 


m 


P 


400 


P 


DO, DO 
D+, D- 
D+ + DO 


(non-leading) 


29, 22 
24, 27 
53 












D- + DO (leading) 
DO + DO + D± 


49 


NA32 


m 


P 


200 


Si 


9 


E743 




P 


800 


P 


DO + DO 


+ D± 


31 


E653 


m 


P 


800 


emulsion 


DO + DO 


+ D± 


146 


WA82 [TTT] 


P 


370 


Si,W 


DO + DO 


+ D± 


266 ± 28 


E789 


gn 


P 


800 


Au 


DO + DO 




~ 2200 


E769 


ins 


P 


250 


Be,Al, Cu,W 


DO + DO 


+ D± + D± 


320 ± 26 


CDF 


P! 


P 


1960 


P 


DO + DO 
D± 

Ds^ 




36 804 ± 409 
28 361 ±294 
5 515 ±85 
851 ± 43 


NA27 


ESI 


n~ 


360 


P 


DO + DO 


+ D± 


57 


NAll 




n~ 


200 


Be 


DO + DO 


+ D± 


29 


NAll 


ED 


7l~ 


200 


Be 


D (leading) 


44 












D^ (non- 


leading) 


30 


NA32 






200 


Si 


D- + DO 
D+ + DO 


(leading) 
(non-leading) 


54 
60 
46 


NA32 






230 


Cu 


DO+DO, D± 


543, 249 
















60, 147 












D- + DO 


+ D*- (leading) 


427 












D+ + DO 


+ D*+ (non-1.) 


425 


E653 


jSH 


7r~ 


600 


emulsion 


DO + DO + D± 


676 


WA75 


m 


7r~ 


350 


emulsion 


DO + DO 


+ D± 


459 


WA82 


[TTT] 


7r~ 


340 


Si,Cu,W 


DO + DO 


+ D± 


2 214 ± 70 


WA82 




7r~ 


340 


Si,Cu,W 


D+ 

D- 




322 ± 20 
449 ± 23 


E769 




7r± 


250 


Be,Al, Cu,W 


DO + DO 


+ D± + D± 


1 665 ± 54 


WA92 


Ei 


7r~ 


350 


Cu,W 


DO + DO 


+ D± 


7 172 ± 108 


E706 


m 


7r~ 


515 


Be,Cu 


D± 




110 


E791 


m 


7r~ 


500 


C 


DO + DO 




88 990 ± 460 



Table 15: Experiments measuring xy and distributions of D mesons. D^ stands for 
D"*" + D~ (and D^ likewise). E789 and CDF only provide px distributions, measured 
in < xf < 0.08 (E789) or in \y\ < 1.0 and px > 5.5 (or 6) GeV/c (CDF). 
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Figure 19: (left) and xp (right) distributions, as measured by the E791 Collabo- 
ration in 500 GeV tt^-C collisions and as calculated by Pythia, version 6.326. 

an even lower PARP(91) value, 0.44 GeV/c. Also the default value of the c quark mass 
has changed, from 1.35 GeV/c^ in version 5.7 to 1.5 GeV/c^ in the latest versions. 
The data points in Fig. IT^ were placed at x values weighted by the function given in 
the original publication jlOj, rather than at the bin centre, following the procedure 
explained in Ref . jll5j . The vertical lines at the top of the figures indicate the bin 
edges. The Pythia curves were scaled up by 2.6, the K-factor extracted from the fit 
to the neutral D meson measurements (see Table ITUI) . 

The highest statistics p\ and xp distributions available for proton induced colli- 
sions were collected by the E769 Collaboration |112j . at E'lab = 250 GeV. The event 
samples are rather small, nevertheless, and the distributions shown in Fig. ISOladd to- 
gether all measured D mesons (D'', D+, and corresponding anti-particles), in data 
sets taken with Be, Al, Cu and W targets. The measured distributions are compared 
with calculations done with Pythia, version 6.326, using the default settings, except 
for PARP(91), which was set to 0.5, 1.0 and 1.5 GeV/c (and PARJ(13) = 0.6). The 
K-factor was set to 3.8 (see Table ITTj) . 

Figure |^ shows the D°, D"*" and p^ differential cross-sections measured by 
CDF in pp collisions at ^/s = 1.96 TeV. Both D and D mesons contribute to these 
single D meson differential cross-sections. The measurements are compared with 
Pythia curves, calculated with three different sets of PDFs, PARP(91) = 1 GeV/c 
and PARJ(13) = 0.6. The K-factors globally fitted to the three data sets are 3.3 
(CTEQ6L), 5.8 (MRST LO) and 5.5 (GRV LO). In the case of the MRST LO PDFs, 
this represents a 50 % increase of the K-factor with respect to the value fitted from 
the fixed-target measurements, 3.8. This observation might indicate that at Teva- 
tron energies the higher-order diagrams missing in our calculation are relatively more 
important than at the lower energies of the fixed-target data. However, the increase 
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E769, E, ^ = 250 GeV, x > 

lab F 

p-A ^ {D°,D°, D*, D', D*, d;.)X 



: MRSTLO{2001) 

Pythia V. 6.32E 

PARP(91) = 0.5GeV/c 

PARP(91) = 1.0 GeV/c 

PARP(91) = I.SGeVfc 



10- 



8 ,10, 
p2 [GeV^/c^] 



50 



E769, E = 250 GeV 

lab 

P-A ^ {□", D°, D*, D', D* 



D}X 



Pythia 6.326 




PARP(91) 


= 0.5 GeV/c 


PARP(91) 


= 1 .0 GeV/c 


PARP(91) 


= 1 .5 GeV/c 




Figure 20: p\ (left) and x-p (right) distributions, as measured by the E769 Col- 
laboration in p-A collisions and as calculated with Pythia, version 6.326, using 
PARP(91) = 0.5, 1.0 and 1.5 GeV/c. 




Figure 21: Single D meson distributions measured by CDF in pp collisions at 
a/s = 1.96 TeV, within \y\ < 1, compared to Pythia calculations. 

is only 10-20 % when we use the CTEQ6L and GRV LO PDF sets, which give pro- 
duction cross-sections with a steeper energy dependence. This exercise shows that 
the K-factors of Table ^2 can be used in the estimates of charm production at the 
intermediate energies of the RHIC experiments, y/s = 200 GeV, where such higher 
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order effects, if any, must be even smaller than at the Tevatron energies. 

In summary, we can say that the existing data on kinematical distributions of 
D mesons can be reasonably well reproduced by calculations made with the latest 
version of Pythia (6.326). 

6.3 D meson pair correlations 

LO diagrams lead to back-to-back cc pairs, with 180° as the difference of azimuthal 
angles of the quarks, A0(cc), and zero pair transverse momentum, pt(cc) |116j . How- 
ever, higher order diagrams, intrinsic and fragmentation lead to smeared A</)(DD) 
and pt(DD) distributions. 

In hadro-production, D meson correlations have been measured in the experi- 
ments NA27 [HI WA75 EHIIIHl, NA32 E653 [1201, WA92 03 [MI and 
E791 |122j . The correlation variables Axp, Ay, A0 and Ap"^ are defined as the 
difference between the D and D values: Axp = xp{D) — xf(D), etc. Also the pair 
variables p|(DD), a;F(DD) and M(DD) have been studied. Out of the quoted ex- 
periments, E791 is the only one which fully reconstructed both D mesons in one of 
the channels D Knir (n=l,2,3). In order to increase the statistical significance of 
their measurements, the other experiments looked for events with secondary vertices, 
characteristic of long-lived particles, irrespective of the reconstruction of the decay 
products. In our comparisons with calculations, we have only used data sets with 
reasonable statistics, as provided by NA32, WA92 and E791. 

The WA92 Collaboration provided A0, Ay, Axp, M(DD), p^(DB) and xf(DD) 
distributions jlHl I121j . for 475 events collected in 7r~-Cu interactions at 350 GeV. 
One of the D mesons, with xp > 0, was fully reconstructed in one of the channels 
D — i> Knvr (n=l,2,3). The second one, which could have any xp, was often only 
partially reconstructed either due to undetected neutral decay products or to the 
limited detector acceptance. To calculate the correlation variables which require the 
momenta of the D mesons, the influence of the neutral decay products was estimated 
by imposing the D meson mass in the reconstruction step and by connecting the sec- 
ondary and primary vertices. A Monte Carlo simulation showed that this estimation 
of the missing information gave a correct calculation of the correlation variables. 

In NA32 both D mesons were only partially reconstructed. A purely topologi- 
cal analysis method selected ~ 500 events with two reconstructed secondary vertices 
well displaced with respect to the primary vertex, resulting in A0, Arj, Ay, M(DD), 
p|(DD) and xf(DD) distributions, for xp > 0, from data collected in 230 GeV 7r~-Cu 
interactions |119j . When there were neutral decays products, an algorithm similar to 
the one used by WA92 estimated the momenta of the D mesons, needed to calculate 
the correlation variables, except for A(f) and At]. The error on the estimated momen- 
tum was evaluated by a Monte Carlo simulation to be ~15%, having little impact 
on the reconstructed correlation variables of the D meson pairs. 

The E791 experiment measured D meson correlations from data collected in 
500 GeV 7r~-C and (a smaller fraction) in 7r~-Pt collisions |122j . reconstructing both 
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D mesons in one of the channels D —>■ Knir (n=l,2,3). The pubhshed A0, Ay, Axp, 
Ap^, M(DD) and p^(DD) distributions are based on 791 ± 44 events. 

We will now compare the three experimental data sets with each other, for each 
of the measured correlation variables. We prepared Pythia curves for each data set, 
simulating 7r~-p collisions at the appropriate energy and properly applying the phase 
space cuts specific of each experiment. The WA92 results were given for full phase 
space, NA32 required that both D mesons were produced at > 0, and E791 
required that both D mesons were produced within —0.5 < y* < 2.5, where y* is the 
rapidity in the centre-of-mass reference frame. 




Figure 22: Acj), pliBB) and distributions measured by NA32, WA92 and E791 
in 7r~ induced collisions, compared to calculations made with Pythia 6.326. 

Figure 1221 shows the transverse correlation variables, A0, p^(DD) and Ap^, mea- 
sured in the three experiments. Since we are only interested in a shape comparison, 
all data sets and respective curves were normalised to each other, taking the E791 
measurement as reference. All measurements are compatible with each other, in 
terms of shape, and can be described reasonably well by the Pythia curves. 

In Fig. 1221 we compare the measured longitudinal correlation variables. Ay and 
xf(DD), with the corresponding calculated curves. These variables are significantly 
affected by the hadronisation step and, therefore, test the modelling of this non- 
perturbative effect. The curves for WA92 (full phase space), depicted as dashed 
lines, are in fair agreement with the measurements. For the E791 case, we show the 
Pythia calculation before (solid line) and after (histogram) applying the phase space 
cuts; the histogram fails to describe the data points. Clearly, the phase space window 
of E791 limits the Ay range to ±3, but the measured distribution seems to be too 
narrow, with respect to Pythia's curve. 

It has been observed by one of the authors of Pythia |123j that the Ay and 
Ax-p distributions measured by E791 could be described quite well (see Fig. |211- 
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Figure 23: Ay and xf(DD) distributions measured by E791 and WA92 in vr" induced 
collisions, compared to calculations made with Pythia for full phase space (solid and 
dashed lines) and after applying the phase space cuts (only for E791, histogram). 




-3 "-2 -1 1 l" d 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Ay >!f "^f 



Figure 24: E791 Ay distribution of DD pairs (left) and WA82 xp distribution of 
D"*" mesons (centre and right) versus Pythia curves calculated with the Peterson 
fragmentation function (left and centre) and with the Lund scheme (right). Figures 
taken from Ref. |123j . 



left) if Pythia's default hadronisation scheme would be replaced by the Peterson 
fragmentation function (with e = 0.05), so that the non-interacting valence quarks 
would not influence the kinematics of the produced charm quark. However, this 
modified model would fail to reproduce the single D meson kinematical distributions 
measured by WA82 jll4j . as shown in Fig. 1241-centre. which are well described by the 
standard scheme (Fig. EH-right). 

In the context of D meson correlations it is worth mentioning the photo-production 
Fermilab experiment "FOCUS" (E831). On the basis of ~ 7000 fully reconstructed D 
meson pairs, they showed |124j that Pythia 6.203 describes very well the correlation 
variables, A0, p^(DD), A|/ and M(DD). Contrary to all other experiments, FOCUS 
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compared the measured data with Pythia events propagated through the simulation, 
reconstruction and analysis algorithms, so that the calculated curves become affected 
by acceptance and efficiencies. This method also takes into account smearing effects, 
and should lead to a more robust comparison between data and theory. In Fig. 
(taken from the FOCUS publication) we can see a remarkable narrowing of the Ay 
distribution, due to the detector acceptance. 
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Figure 25: Correlations for fully reconstructed DD pairs from the photo-production 
experiment FOCUS. For comparison, events generated with Pythia 6.203 (dotted 
lines) were tracked and reconstructed by the simulation code (histograms), being 
affected by acceptances, efficiencies and smearing effects, as the data. Figures taken 
from Ref. ^2M- 




Figure 26: DD invariant mass distributions measured by NA32, WA92 and E791, 
in 7r~ induced collisions, compared to calculations made with Pythia, version 6.326. 
The lines show Pythia calculations prepared for full phase space while the histograms 
show the Pythia curves with the correct phase space cuts applied. 

NA32, WA92 and E791 also measured the DD invariant mass distributions, as re- 
ported in Fig. |2ni While the NA32 and WA92 distributions are nicely reproduced by 
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Pythia's curves, the E791 points show a spectral shape quite different from the calcu- 
lation, with a very unexpected rise when approaching the DD production threshold, 
in apparent conflict with phase space considerations. 

To summarise, the latest version of Pythia (6.326), with default settings (except 
for PARP(91), which should be kept between 1 and 1.5 GeV/c), is able to repro- 
duce reasonably well the available D meson pair correlation measurements, at least 
as provided by the NA32, WA92 and FOCUS Collaborations. However, there is a 
significant disagreement in what concerns the Ay, Axp and M(DD) distributions 
measured by E791, which might be due to a problem in the 8-dimensional acceptance 
correction procedure. The highest statistics data set available for D meson corre- 
lations, collected by FOCUS, is in good agreement with Pythia events propagated 
through the simulation of the apparatus. 

We finish this section by mentioning a different kind of pair correlation, the relative 
abundance of D~ with respect to D"*", evaluated through the asymmetry 

^^^^^-a(D-) + a(D+) " 

A{xf) has been measured in vr" induced collisions by the WA82 |114j . WA92 [HI [121], 
E769 [Sni and E791 |125j Collaborations. In Fig. EH we show a comparison (made 
by the authors of Pythia |lU5j ) between the available data and a Pythia calculation 
specifically made for 340 GeV 7r~p collisions. The solid line, labelled "pair produc- 
tion", corresponds to a calculation only including the standard LO processes of charm 
production, plus initial and final state radiation, while the dashed line labelled "all 
channels" also includes "flavour excitation" and "gluon splitting" diagrams. 




Figure 27: Asymmetry of charm particle production. Figure taken from Ref. |l()5j . 

These measurements, and the calculations, clearly indicate that the D~ mesons 
are produced at more forward rapidities than the D+ mesons, confirming that the D~ 
is the leading particle in vr" induced collisions. 
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6.4 Nuclear dependence of charm production 

Most experiments used nuclear targets and published the D meson cross-sections for 
p-N and vr-N collisions, assuming a linear dependence with the mass number of the tar- 
get nucleus, A. However, as explained in Section we expect that ant i- shadowing 
in the parton distribution functions will affect charm production at energies below 
^/s ~ 150 GeV. Therefore, the p-N and vr-N values obtained by extrapolating the p-A 
and TT-A measurements using the linear A scaling should be somewhat higher than 
the values directly measured in pp and vr-p collisions. From Fig. |H1 we can see that 
the nuclear anti-shadowing given by the EKS98 model increases with energy between 
^/s = 20 and 40 GeV, the energy range where the data have been collected. But even 
at ^/s = 40 GeV, and for the Pb nucleus, the expected increase in the charm produc- 
tion cross-section due to the nuclear anti-shadowing effect, with respect to the linear 
extrapolation of pp collisions, is only around 10%, a value too small to be visible in 
the presently available measurements, given their rather large uncertainties. 



Exp. 


Eiab[GeV] Target 
Phase space 


Observed 
D mesons 


a 


p-A collisions 


E789 

m 


800 Be,Au 
0<xf<0.08, pt<1-1 


Be: 1360 
Au: 1040 DO 


DO : 1.02 ±0.03 ±0.02 


TT -A collisions 


WA82 

m 


340 Si,Cu,W 
xf > 0.0 


Si: 102 (D0,D+) 
Cu: 528 (DO, D+) 
W:1017 (D0,D+) 


DO+D+ : 0.92 ±0.06 
DO ^ Kvr : 1.03 ±0.11 
DO ^ Kmrn : 0.93 ±0.11 
D+ Knn : 0.84 ± 0.08 


E769 

m 


250 Be,Al,Cu,W 
Xf > 0.0 


all targets: 
650 DO 
776 D+ 


DO+D+ : 1.00 ±0.05 ±0.02 
DO : 1.05 ±0.15 ±0.02 
D+ : 0.95 ±0.06 ±0.02 


WA92 

m 


350 Cu, W 
Xf > 0.0 


Cu: W: 
3245 DO, 628 DO 
2753 D+, 546 D+ 


DO+D+ : 0.93 ±0.05 ±0.03 
DO : 0.92 ±0.07 ±0.02 
D+ : 0.95 ±0.07 ±0.03 


E706 

m 


515 Be, Cu 
xf>-0.2, 1<pt<8 


Be+Cu: 110 D+ 


D+ : 1.28 ±0.33 



Table 16: Nuclear target dependence in proton and pion induced collisions. Note 
that DO and D"*" mean DO + DO and D"^ ± D~, respectively, in GeV/c. 



Some experiments made measurements with two or more targets and fitted their 
data to the A'^ form, extracting the a values compiled in Table El and Fig. |2H1 In 
the third column we give the number of observed D mesons for each target, except if 
only the total number was provided. The errors of the WA82 measurements include 
systematic uncertainties, which are small with respect to the statistical error. All the 
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• neutral D 
□ charged D 


E789 




2400 


Be 


Au 


WA82 




1647 


Si, Cu, W 


E769 




— 650 
776 


Be, Al, Cu, W 


WA92 


— 1 


3873 
3299 


Cu, W 


E706 






Be 


Cu 






0.5 


1 


1.5 




2 a 



Figure 28: Values of the a parameter extracted from data collected in pion and proton 
induced collisions. The two vertical lines indicate a ±10 % range around a = 1. 

published values are consistent with a = 1, within their large errors. 

In a dedicated study |126j , using the relative cross-sections of four nuclear targets 
(Be, Al, Cu and W), E769 extracted a = 1.00 ± 0.05 ± 0.02 as the best description 
of the measurements, made with a pion beam. Within the statistical accuracy of the 
data, they saw no change between the values obtained with the 7r~^ and 7i~ beams, or 
any dependence of a on the D meson's pi and xp. Also WA82 and WA92 have not seen 
any dependence of a on xp or p^, but with even larger uncertainties. The existence 
of such a dependence (seen for many other particles, including the cannot be 
excluded from the present measurements and, if observed by a future high statistics 
experiment, could imply a revised extrapolation of the E706 and E789 measurements 
to elementary vr-N/p-N cross-sections, given their limited pt coverages. 

The WA78 Collaboration also studied the nuclear dependence of the charm pro- 
duction cross-section, by measuring the yield of prompt single muons in a beam- 
dump experiment, both with 320 GeV vr" and with 300 GeV proton beams, in- 
teracting on Al, Fe and U targets. The resulting a values (for xp above ~0.1) 
were a(/i+) = 0.79 ± 0.12 and a(/i-) = 0.76 ± 0.13 for the p-A data [HTl, and 
a(/i+) = 0.76 ± 0.08 and a(/i-) = 0.83 ± 0.06 for the tt-A data [HHl. These results 
are significantly lower than 1, raising doubts on whether the prompt single muons 
can be cleanly ascribed to semi-leptonic decays of charmed particles. 

We close this discussion with a final remark concerning the cXpA = ctq ■ param- 
eterisation. Figure El collects the D meson production cross-sections measured in 
proton and pion induced collisions. Three measurements were made with hydrogen 
targets (NA16, NA27 and E743) while the other values were measured with nuclear 
targets and divided by A. Considering the error bars of the data points, we do 
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not see any significant difference between the a^p^j A values and the cXpp values. In 
Ref. |129j (from 1987) it was argued that the elementary charm production cross- 
sections derived from p-A measurements (performed with nuclear targets of A > 9) 
should be parameterised as UpA = -ft'o ■ c^pp ■ following observations made with 
light flavour data |13Uj (where Kq is around 1.5-2.0). Is was noted, in particular, 
that = 1.5 would result in a = 1.0 at xp = when comparing NAll (vr-Be) and 
NA27 (vr-p) data |129j . Hence, the dpA/A values should be higher than the pp values, 
or a^'^ 7^ (Tpp, as recently recalled in Ref. |131j . However, looking at the error bars in 



we see that Kq = 1.5 implies a = 1.0!I^q4, ^ very poor statistical 

1.2ln^, perfectly compatible 



-0.4) 



Fig. 1 b of Ref. 

significance. In fact, setting Kq = 1.0 would lead to a 
with 1.0, within errors. . . The claim of Ref. |129j is, thus, unsubstantiated (especially 
when we note that systematic uncertainties were neglected in the comparison between 
the data of the two experiments). 
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Figure 29: Quality of the fit to the charm production cross-section data shown in 
Fig. El as a function of the Kq factor used to scale the pp points. 

We have revisited this issue using the charm production data currently available, 
almost 20 years later. Would the new measurements prefer a Kq value significantly 
higher than 1? We fitted the charm production cross-sections shown in Fig. IT^ after 
scahng the pp values by a Kq factor. Figure EHl shows how the fit quality, expressed in 
terms of per degrees of freedom, changes when varying Kq. We see that the best 
fit is obtained with Kq ~ 1.1, and that the fit quality significantly degrades when 
using Kq ~ 1.5. 



6.5 Beauty production cross-sections 

As we have seen in Sectional not many experiments have measured beauty production 
cross-sections, in pion or proton induced collisions, and most of those measurements 
are derived under mo del- dependent assumptions. This might partially explain why 
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some data points, collected at essentially the same energy, differ by factors of 5 (NAlO 
and WA78; E771 and E789). The available measurements are shown in Fig. 1201 as a 
function of ^/s, separated between pion (left) and proton (right) data. Two points are 
shown for NAlO, using open and closed squares, corresponding to different theoretical 
assumptions for the kinematical distributions used to derive the full phase space value. 




Figure 30: BB production cross-sections in pion (left) and proton (right) collisions. 

These measurements are compared to the results of the calculations we have made 
with Pythia. On the left panel we show the curves calculated with GRV-P LO (pion) 
and MRST LO (proton) PDFs, scaled with K-factors of 1, 2 and 3. Given their 
significant spread, it is not meaningful to use the data points to fit a specific K-factor 
value. On the right panel we show the fixed-target proton data points, together with 
pp data from UAl and CDF. The two curves were calculated with CTEQ6L and 
MRST LO PDFs, and were (arbitrarily) scaled up with a K-factor of 2. We remark 
that the beauty production cross-section at such high energies is dominated by gluon 
fusion and, therefore, is identical for pp and pp collisions. 

It is clear that the fixed-target points, on their own, are not able to give a mean- 
ingful normalisation of the calculations, and we must profit as much as possible from 
the three measurements provided by the pp collider experiments: UAl quotes the 
full phase space cross-section at = 630 GeV, extrapolated from > 6 GeV/c 
and \y\ < 1.5, while CDF gives cross-sections at ^/s = 1.8 and 1.96 TeV, within their 
specific phase space window. 

In Run I, CDF measured a single B"*" meson cross-section of 3.6 ± 0.4 ± 0.4 fih 
(dividing the sum of reconstructed B"*" and B~ mesons by 2), in the kinematical 
window \y\ < 1.0 and > 6.0 GeV/c. According to Pythia (version 6.326), this 
window covers 7.8 % of the full phase space. The corresponding calculated cross- 
section is 1.9 /ib if evaluated with the CTEQ6L PDF set, and 1.2 /ib if evaluated with 
the MRST LO or GRV LO PDF sets. In Run II, CDF determined the single B hadron 
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(beauty only; not anti-beauty) production cross-section in \y\ < 0.6, i7.6iO.4l2J A^b, 
using displaced J/V^'s tagged to come from beauty hadron decays. In this kinematical 
window (~22% of the full phase space), Pythia gives 16, 9 and 10 /ib (corresponding 
to K-factors of 1.1, 1.9 and 1.7), with the CTEQ6L, MRST LO and GRV LO PDF 
sets, respectively. These two measurements are also shown in Fig. EOl and should be 
compared to the Pythia curves calculated for the corresponding phase space windows. 

Taking into account that the E771 and E789 data points (at the same energy) 
differ by a significant amount, and that the collider values have been measured in 
kinematical windows covering only a relatively small fraction of full phase space, it 
is remarkable that the Pythia curves, with a K-factor of 2, go through essentially all 
the data points within around a factor 2, over four orders of magnitude in beauty 
production cross-section, between the fixed-target and the CDF energies. It is partic- 
ularly remarkable that the comparison between the two most reliable measurements, 
HERA-B and CDF, does not indicate any increase of the K-factor with energy. 

Using the Pythia curves we derive a beauty production cross-section in pp col- 
lisions at RHIC energies, ^/s = 200 GeV, of around 2.5 /ib. This value is in good 
agreement with the prediction of a QCD FONLL calculation, which goes beyond the 
("fixed-order") NLO result by including the resummation of next-to-leading loga- 
rithms ("NLL"). See Ref. jHE] and references therein for further details on the calcu- 
lation and on its "parameters": the heavy quark mass, the strong coupling, Aqcd, the 
PDF set, the factorisation and renormalisation scales, the fragmentation functions, 
etc. The result of this calculation, a{hh) = 1.87to;g7 /^b, has been included in Fig.lHIH 
for comparison purposes. The error bar represents the uncertainty of the calculation, 
estimated by varying certain input parameters within reasonable ranges. 

6.6 Beauty kinematical distributions 

Measurements of kinematical distributions of beauty mesons are rare. UAl pub- 
lished a pt distribution on the basis of around 3000 events (merging four indepen- 
dent analyses), for > Q GeV/c and \y\ < 1.5 jH]. From Run I, CDF published 
a pt distribution for B+ mesons with p^ > Q GeV/c and \y\ < 1.0, based on 387 
events Using around 40 pb~^ of Run II data, CDF tagged ~ 40 000 J/ip mesons, 
within \y{J/ip)\ < 0.6, as being produced away from the interaction vertex, presum- 
ably by beauty decays The corresponding p^ distribution is shown in Fig. 1311 
where it is compared with Pythia calculations made with three PDF sets and with 
PARP(91) = 1 GeV/c, on the left panel, and with a QCD FONLL calculation [HSI, 
on the right panel. 

The agreement between the calculations and the measured data is quite good, 
over three orders of magnitude and down to very low values, where theoretical 
uncertainties are particularly important. The curves calculated with Pythia were 
normalised using the K-factors derived in the previous section. 
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Figure 31: pt distribution of J/ip's from beauty decays, as measured by CDF, within 
\y\ < 0.6. The error bars include an overall 6.9% systematic uncertainty |132j . The 
curves on the left panel were calculated with Pythia, version 6.326; those on the right 
panel are a FONLL calculation |133j (figure taken from Ref. |55j). 

6.7 Beauty feed- down to J/?/; 

The study of the decay mode B — * J/ipX allowed several experiments to determine 
the beauty production cross-section. On the other hand, for studies of J/ip production 
this decay channel constitutes a source of background which should be carefully eval- 
uated. This feed-down source of J/ip mesons is particularly important in the context 
of the study of J/ip suppression (or enhancement) in heavy-ion collisions at suffi- 
ciently high energies, as those at the RHIC and LHC colliders. Since beauty mesons 
should not be affected by the medium formed in heavy-ion collisions, the fraction of 
J/ip^s which come from beauty decays will not be suppressed. Therefore, the pattern 
of J/ip suppression measured at RHIC or LHC energies cannot be directly compared 
with predictions based on J/ip melting by QGP formation, for instance, without con- 
sidering the beauty "contamination" . This observation emphasises the importance of 
measuring the beauty yield in the RHIC and LHC heavy-ion experiments. 

We will now make a rough evaluation of the level of this problem. Figure IH^ shows 
the fraction of J/ip^s coming from the decay of beauty hadrons, as a function of the 
J/V^'s transverse momentum, as measured by CDF in Run II, within \y\ < 0.6. We 
see that the fraction of J/ip^s from beauty decays rises steeply with increasing px, 
from ~ 10 % at < 3 GeV/c to ~ 50 % at pt ~ 20 GeV/c. At = 630 GeV, UAl 
observed that 31 ± 2 ± 12% of all J/^/^'s produced within 5 < ^ 25 GeV/c and 
\y\ < 2.0 come from beauty decays jSHl- This fraction becomes negligible (less than 
0.1 %) at the much lower energy of HERA-B [FTj . 

Since beauty production is expected to scale linearly with the mass number of 
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Figure 32: Fraction of J/?/) mesons from beauty decays as a function of the J/^/^'s 
transverse momentum, as measured by the CDF experiment for \y{^/%l))\ < 0.6, at 
^/s = 1.96 TeV jSSj- Error bars include statistical and systematic uncertainties. 



the colliding nuclei, while J/ip production scales less than linearly (as at SPS 
energies |134p . in Au-Au or Pb-Pb collisions the relative fraction oi J/ip mesons 
resulting from beauty decays will be higher than in pp collisions. If direct J/ip 
production is further suppressed in heavy-ion collisions (NA50 measured a factor 2 
of extra suppression in central Pb-Pb collisions at the SPS |135p . beauty production 
could account for a significant fraction of the observed J/ip yield, especially at LHC 
energies and at high values. Studies of J/ip suppression as a function of pt are 
particularly sensitive to this feed-down source, given its strong p^ dependence. 



6.8 Charm cross-section measurements without vertexing 

Besides the charm production cross-section measurements mentioned in Section 0J 
mostly performed with especially designed detectors and affected by relatively low 
background levels, there are a few measurements made by other experiments, in 
more difficult conditions. We will consider in this section three of these "indirect 
measurements" , recently made by experiments working on the field of "quark matter 
physics": NA50 at the SPS; PHENIX and STAR at RHIC. While these experiments 
are mostly devoted to the study of high-energy nuclear collisions, they have also taken 
pp, p-A or d-Au data. We will only consider results from these more elementary 
collisions. 

The NA38 and NA50 Collaborations studied dimuon production, of mass above 
1.5 GeV/c^, in p-A, S-U and Pb-Pb collisions, at SPS energies |136j . to look for 
evidence of thermal dimuon production from a quark-gluon plasma, presumably 
formed in heavy-ion collisions. The measured mass distribution, in the continuum 
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surrounding the ^/ip and if)' resonances, was compared to the superposition of two 
expected sources: dimuons from the Drell-Yan mechanism and muon pairs from si- 
multaneous semi-muonic decays of pairs of D mesons. The proton-nucleus data, 
collected with a 450 GeV proton beam and several nuclear targets (Al, Cu, Ag, W), 
could be rather well described, both in terms of mass and distributions, as the 
sum of these two contributions, simulated with Pythia 5.7, with MRS A PDFs and 
PARP(91) = 0.9 GeV/c. Pythia was also used to calculate the "extrapolation factor" 
needed to go from the elementary (nucleon-nucleon) charm production cross-section 
used in the event generation to the yield of muon pairs from D meson pair decays 
detected in the phase space window probed by the experiment: 3 < yiab < 4 and 
I cos(^cs)| < 0.5, where yiab is the dimuon rapidity in the laboratory frame and ^cs 
is the CoUins-Soper polar angle |137j . Using this calculated factor and a global fit to 
the four p-A dimuon mass distributions, NA50 derived a'^f = 36.2 ±9.1 /ib as the 
full phase space charm production cross-section, per nucleon, which best reproduces 
the dimuon data collected in p-A collisions at 450 GeV. The derivation used a lu- 
minosity deduced from the number of J/ip events (and a previously measured J/ip 
cross-section), and assumed a linear dependence of a^c with the target mass number. 




Figure 33: Left: The dimuon mass and cos(^cs) variables are strongly correlated for 
muon pairs from charm decays. Right: Dimuon mass distributions for Drell-Yan and 
charm decays, before (thick lines) and after (thin lines) applying the | cos(6'cs)| < 0.5 
cut. 

We should be cautious when comparing this value with those directly obtained 
from the study of hadronic decays of D mesons. The "charm yield" is extracted from 
the dimuon mass spectra after subtracting a very important "combinatorial back- 
ground" of muon pairs due to pion and kaon decays, and a Drell-Yan contribution 
which, at such low dimuon masses (1.5-2.5 GeV/c^), is rather uncertain. Further- 
more, the relative contribution of open charm decays and Drell-Yan dimuons to the 
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mass distribution is quite sensitive to the kinematical cuts imposed on the events, 
given the strong correlation between the dimuon mass and the cos(^cs) kinematical 
variables. Figure l^^left shows that when we select high mass muon pairs from DD 
decays we are forcing them to have cos(6'cs) close to —1 or +1. Such events can- 
not be measured by the NA38/50 spectrometer, which imposes the kinematical cut 
I cos(^cs)| < 0.5. The influence of this cut on the mass distributions of muon pairs 
resulting from charm decays and on Drell-Yan dimuons is shown in Fig. I^^right. 
for the case of pp collisions at -Eiab = 400 GeV. We see that the Drell-Yan dimuon 
mass distribution is simply scaled down while the reduction of the charm spectrum 
is much more pronounced for the higher masses, leading to a significant change of 
the shape. It is curious to note, in particular, that the Drell-Yan contribution would 
not dominate the mass range above the J/V' peak (a common assumption at these 
collision energies) if the | cos(^cs)| < 0.5 cut would not be applied. 

The PHENIX Collaboration measured inclusive single electron spectra in pp 
collisions ai = 200 GeV [HHI, using data collected in RHIC Run-2 (2001/2002), 
in the following phase space window: 0.4 < < 5.0 GeV/c, |?7| < 0.35, A0 = 7r/2. 
Charged particles were tracked in a drift chamber and a pad chamber layer. Electrons 
were selected using the information provided by an electromagnetic calorimeter and 
a ring imaging Cerenkov detector. The ratio of "non-photonic" electrons, assumed to 
come from charm and beauty decays, to background (photonic) sources is ~0.4 for 
Pt < 1-5 GeV/c and > 1 for higher values. The shape of the spectrum obtained 
after background subtraction is reasonably well described by a superposition of charm 
and beauty contributions, as simulated with Pythia, in the range up to 1.5 GeV/c. 
In the higher px range the measured distribution exceeds the calculation, made with 
Pythia 6.205, with = 1.25 GeV/c^ s as the definition and the CTEQ5L 
set of PDFs, as described in Ref. |110j . Since the total production cross-section is 
dominated by the low px region, the high-px disagreement can be neglected, and the 
normalised Pythia charm curve was integrated down to px = 0. The resulting mid- 
rapidity cc production cross-section is dacc/dy = 0.20 ± 0.03(stat) ± O.ll(syst) mb 
and the estimated full phase space value is Ucc = 0.92 ± 0.15 ± 0.54 mb. 

The STAR Collaboration extracted the cc production cross-section in d-Au colli- 
sions at ^/s = 200 GeV from data collected in 2003 jl39j . Two independent measure- 
ments were made: using inclusive electron data to probe charm semi-leptonic decays 
and using directly reconstructed neutral D mesons, through the D —* Kir hadronic 
decay channel. The yield of neutral D mesons with \y\ < 1 and px in the range 
0.1-3 GeV/c was extracted from the invariant mass spectrum of kaon-pion pairs, 
reconstructed in the TPC. The lack of accurate vertexing information imposed the 
pairing of all oppositely charged kaon and pion tracks of each event, resulting in a 
signal-to-background ratio around 1/600. After combinatorial background subtrac- 
tion, by event mixing, the D° yield was determined with an uncertainty estimated to 
be around 15%. Electrons and positrons were identified combining the dE/dx mea- 
sured in the TPC with the velocity information provided by a TOF system (covering 
~ 7r/30 and —1 < r] < 0). After subtracting the estimated contribution from 
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photonic sources, the electron distribution in the range 1-4 GeV/c was compared 
to calculations of the yields expected from charm decays. The D'' yield obtained 
from a combined fit to the D'' and electron data was converted in a full phase space 
charm production cross-section, using a pp inelastic cross-section of cr?*^) = 42 mb, 
the average number of binary nucleon-nucleon collisions in d-Au collisions (7.5), a 
factor of 4.7 ±0.7 to extrapolate from mid-rapidity to full phase space, and a fraction 
of D° mesons with respect to all charmed hadrons of 0.54 ± 0.05 (measured at LEP). 
The result was 1.4 ±0.2 ±0.4 mb per nucleon-nucleon interaction. While this value is 
considerably higher than the PHENIX pp cross-section, obtained at the same energy, 
the two results are compatible with each other within (systematic) uncertainties. 

In Fig. IH^left we compare the three measurements mentioned above, made by 
NA50, PHENIX and STAR, to the calculations we have made with Pythia and already 
presented in Fig. ^| The curves were normalised with the K-factors determined from 
the fixed-target measurements, and reported in Table ^1 Taking as reference the 
calculations made with the CTEQ6L set of PDFs, the value derived from the NA50 
p-A dimuon data is a factor 2.1 ± 0.5 too high. The PHENIX and STAR values are 
closer to the calculations, differing by factors of 1.15 ± 0.7 (PHENIX) and 1.7 ± 0.6 
(STAR). With respect to the calculation made with MRST LO PDFs, the NA50 
value remains a factor 2 too high while the PHENIX and STAR points become too 
high by factors 2.1 ± 1.3 (PHENIX) and 3.2 ± 1.0 (STAR). 
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Figure 34: Left: Comparison of "indirect" (Tcc measurements with Pythia calculations. 
Right: Influence of the rric, and PDFs settings used by PHENIX on the calculated 
cc production cross-sections. 

We have also included in this figure, for comparison with the RHIC data points and 
with Pythia's curves, the result of the FONLL calculation for charm production 
at ^/s = 200 GeV: cr(cc) = 2561^4° A^t). The rather large uncertainty was evaluated 
by varying the input settings of the calculation (see the text describing the FONLL 
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point of Fig. 1201 for further details). 

On the right panel of Fig. we show four different curves, all calculated with 
Pythia 6.326 but with different settings. The solid line is the curve obtained with 
our standard settings and CTEQ6L PDFs, already presented on the left panel. The 
dashed line is obtained with the settings used by PHENIX: rric = 1.25 GeV/c^, s 
as the definition and the CTEQ5L set of PDFs. The dotted line uses "interme- 
diate" settings: the same mass and definition as used by Phenix, but using the 
more recent CTEQ6L PDFs. This allows us to visualise the differences between the 
settings in a more gradual way. Finally, the dashed-dotted curve is done with the 
same "standard" settings as the solid line, except that we used CTEQ6M instead of 
CTEQ6L PDFs. We recall once more that the use of NLO PDFs, like CTEQ6M, is 
not appropriate in LO calculations, such as those performed by Pythia. 

It is a curious coincidence that the solid and dashed curves cross each other at 
the energy of the RHIC data points, in spite of the fact that only the fixed-target 
measurements contribute to their normalisations. The change in c quark mass, in 
definition and in PDF sets compensate each other, so that both the default settings, 
with CTEQ6L, and the "Phenix settings" give a cc cross-section at y/s = 200 GeV 
of 800 yub. The difference between the dashed and dotted lines is due to changes 
between the CTEQ5L and CTEQ6L PDFs, and probably originate in the harder 
gluon densities of CTEQ6L: less gluons at low x (see Fig. E)) results in somewhat 
lower cross-sections at high energies. 

All the calculations give essentially the same values in the SPS energy range, 
since they are "forced" by the existing measurements to cross each other at around 
^/s = 30 GeV (see Table IT^ . In the context of the NA50 and NA60 experiments, 
however, it is worth recalling that at these energies charm production is expected to be 
sensitive to nuclear effects on the PDFs, according to the EKS 98 model (see Fig. |H)). 
For i?iab = 400 GeV, the charm production cross-section per nucleon is expected to 
be enhanced by ~ 4 % in p-Be and by ~ 10 % in p-Pb collisions. In the case of Pb- 
Pb collisions at E\sh = 158 GeV per nucleon, the corresponding enhancement factor 
should be around 12 %. 

7 Summary and conclusions 

In this report we reviewed some aspects of open charm and beauty hadro-production 
from data collected in fixed-target and collider experiments. We considered measure- 
ments made with proton and pion beams, at the CERN SPS, DESY and Fermilab, 
or in pp colliders, at CERN (UAl) and Fermilab (CDF). 

We used the Monte Carlo event generator Pythia, version 6.326, as a theoretical 
model to describe the reviewed data. We have seen that Pythia provides a rea- 
sonable description of the available charm data, in terms of energy dependence of 
production cross-sections, D meson kinematical distributions, and pair correlations, 
provided the intrinsic transverse momentum of the colliding partons, k^, is generated 
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with a Gaussian of width PARP(91) ~ 1 GeV/c. However, we observed that the ra- 
tio between charged and neutral D meson production cross-sections, as measured in 
proton induced colhsions, is around a factor of 2 higher than the value expected from 
Pythia's calculations, with default settings. The comparison with data collected in 
other collision systems, and at many different energies, shows that naive spin count- 
ing arguments do not seem to apply in the case of D meson production, presumably 
because of the significant mass difference between the D and the D* states. This 
indicates that Py (parameter PARJ(13) in Pythia) should be set to 0.6 (when gen- 
erating charm events; for beauty the default value of 0.75 should be kept). Except 
for PARP(91) and PARJ(13), we used Pythia (6.326) with default settings. Besides, 
the calculations need to be up-scaled by an appropriate K-factor. In particular, using 
nic = 1.5 GeV/c^ and the CTEQ6L set of PDFs, the best description of the measured 
charm cross-sections requires a K-factor of 3. 

Normalising the calculated charm production cross-sections to the existing (fixed- 
target) data, we have deduced the values expected for the RHIC energies, and com- 
pared them with the values extracted from the measurements of PHENIX and STAR. 
The issue of nuclear effects in the parton densities was also addressed in some detail. 
At SPS energies, the EKS 98 parameterisation leads to a ~10% higher cc cross- 
section in p-Pb collisions, with respect to a linear extrapolation from proton-proton 
collisions. Such a small effect cannot be confirmed by looking at existing data, given 
the very large uncertainties of the measurements. Maybe further insight into these 
issues will result from the study of the proton-nucleus data collected by the NA60 
experiment in 2004, with 400 GeV protons incident on seven different nuclear targets 
simultaneously placed on the beam line. 

There are not many measurements of beauty production cross-sections and most 
of those obtained with pion beams are significantly "model dependent" (especially 
the oldest values). In these conditions, it is remarkable (even surprising) that the cal- 
culated energy dependence of the beauty cross-section, for pp or pp collisions, is able 
to reproduce rather well the HERA-B and CDF measurements, which differ by four 
orders of magnitude. This results in a relatively accurate prediction for the beauty 
production cross-section in the "intermediate energies" of the RHIC experiments, 
2.5 /ib. Nevertheless, we must emphasise the relevance of a direct measurement of 
the beauty production cross-section at RHIC, particularly important to ensure a cor- 
rect interpretation of the J/ip measurements, in heavy- ion collisions, as a function of 
collision centrality and of pt- 
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Note added in proof 

After this paper was completed, the HERA-B Collaboration finished their open charm 
analysis |14()j . The resulting production cross-sections, per target nucleon and in the 
phase space window —0.15 < xf < 0.05, are 10.7 ± 1.2 ± 1.4 fih for the charged D 
mesons and 26.3±2.4±2.6 fih for the neutral ones. These values must be extrapolated 
to full phase space before they can be compared to other measurements. According 
to Pythia, version 6.326, the HERA-B xp window covers slightly different fractions of 
full phase space for the leading and non-leading D mesons, because they have different 
rapidity distributions: 49 % for the D"*" and D°; 58 % for the D~ and D^. These values 
were calculated with the CTEQ6L set of PDFs and would be somewhat smaller if 
calculated with the MRST LO (or GRV LO) PDFs: 43% and 55%. Averaging the 
extrapolation factors obtained with these three PDF sets and taking into account 
the different particle and antiparticle values, the total production cross-sections are 
21.4 ± 2.4 ± 2.8 yub and 52.6 ± 4.8 ± 5.2 /xb for the charged and neutral D mesons, 
respectively. These values replace the preliminary ones given in Table |3 




Figure 35: Neutral (closed symbols) and charged (open symbols) D meson cross- 
sections (left) and corresponding total cc cross-sections (right), as a function of y/s, 
compared to Pythia calculations made with three different PDF sets. 
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The left panel of Fig. 1351 shows the neutral and charged D meson total production 
cross-sections, as a function of a/s, including the final HERA-B values, compared to 
curves calculated with Pythia with Py = 0.6 and three different LO PDF sets (this 
figure supersedes the top panels of Fig. [T21). The right panel of this figure updates 
the total cc cross-sections, previously shown in Fig. ITBl-left. replacing the preliminary 
HERA-B value (given in TablelHI) by the new value: 44.4±3.2±3.5 /xb. The K-factors 
resulting from the new fits to the data points are 2.8, 3.6 and 4.4 for the CTEQ6L, 
MRST LO and GRV LO PDF sets, respectively (these values supersede those given 
in Table ini). 

It is particularly remarkable that the final HERA-B value for the ratio of charged 
to neutral D meson cross-sections is 0.40 ± 0.06 ± 0.04, in very good agreement with 
the value expected when setting Py = 0.6. Figure EHl shows the updated version of 
the left panel of Fig. [T21 
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Figure 36: Ratio between the charged and neutral D meson cross-sections, for proton 
induced collisions, updated with the final HERA-B value. 
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